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SUMMARY OF METHODS FOR CALCULATING DYNAMIC LATERAL STABILITY AND RESPONSE 
AND FOR ESTIMATING LATERAL STABILITY DERIVATIVES ! 


Ву донм P. CauPBELL and Marion О. McKiINNEY 


SUMMARY 


А xummary of methods for making dynamic lateral stability 
and response calculations and for estimating the aerodynamic 
stability derivatives required for use in these calculations is pre- 
sented. The processes of performing calculations of the time 
histories of lateral motions, of the period and damping of these 
motions, and of the lateral stability boundaries are presented as 
a series of simple straightforward steps. Existing methods for 
estimating the stability derivatives are summarized and, in 
sume cases, sumple new empirical formulas are presented. Ref- 
erence is also made to reports presenting experimental data 
that should be useful in making estimates of the derivatires. 
Detailed estimation methods are presented for low-subsonic- 
speed. conditione but only а brief discussion and a list of refer- 
ences are given for transonic- and supersonic-apeed conditions. 


ENTRODUCTION 


Dynamic lateral stability has not received widespread 
attention in the past because it has not generally been a 
serious problem in the design of airplanes. Consideration 
of dynamic lateral stability has recently become more im- 
portant, however, because current design trends toward the 
use of low aspect ratio, sweepback, and higher wing load- 
ing have, in many сазев, led to unsatisfactory dynamic lat- 
eral stability. Airplane designers are therefore finding 1t 
necessary to make such calculations in connection with the 
deaign &nd modifieation of airplanes. In many cases these 
calculations are difficult to perform for designers who have 
had no previous experience in theoretical stability work 
because most of the published theoretical analyses are not 
presented in a form that is especially suited to Ше compu- 
tation of dynamic stability. The estimation of the stability 
derivatives required in dynamic stability calculations has 
also been found to be difficult in many cases. Although 
theoretical and experimental data on these derivatives have 
appeared in numerous publications, no single publication has 
presented methods for estimating the derivatives for all 
types of airplanes. 

One approach to а presentation of methods of calculating 
stability and estimating stability derivatives in a form suit- 
able for use by designers was made by Zimmerman in 
reference 1. Although this report has proved to be of valu- 
able assistance to designers in making dynamic stability 
calculations, recent trends in airplane design have caused its 


usefulness to be seriously limited. For example, the equa- 
tions of reference 1 do not include Ше product-of-inertia 
terms which have been shown by recent studies to be very 
important in some cases. (See references 2 and 3.) More- 
over, the calculation of the time histories of lateral motions, 
one type of calculation that has been the subject of increasing 
interest in the last few years (references 4 to 7), is not covered 
in reference 1. The methods of estimating stability deriva- 
tives presented in reference 1 are also limited because they 
apply only to airplanes having unswept wings with an aspect 
ratio of 6 and operating at speeds at which compressibility 
effects are negligible. The purpose of the present report is 
to extend the methods of reference 1 to include the methods 
of computation which are of current interest to designers 
and to include methods of estimating derivatives for con- 
figurations and flight conditions which are now being 
considered. 

This report summarizes and reduces to simple straight- 
forward steps methods for computing the time histories of 
lateral motions, the period and damping of these motions, 
and the lateral stability boundaries. Existing methods of 
estimating stability derivatives for a variety of airplane 


configurations are summarized and, in some cases, simple — 


new empirical formulas are presented. Reference is also 
made to reports presenting experimental data that should be 
useful in making estimates of these derivatives. 


SYMBOLS 


All forces and moments are referred to the stability 
system of axes which is defined in figure 1. The following 
definitions apply to the symbols except where they are 
otherwise defined: 

т mass of airplane, slugs 

© wing area, square feet 

c wing mean chord, feet (5/2) 
b wing span, feet 


Va span of that part of wing that has tip 
dihedral, feet 
l tail length (distance from center of 


pressure of vertical tail to center of 
gravity, measured parallel to longi- 
tudinal stability axis; values of l 
must be calculated for each angle 
of &ttack), feet 


t Supersedes NACA TN 2400, “Summary of Methods for Calculating D ynamic Lateral Stability and Response and for Estimating Lateral Stability Derivatives” by John P. Campbell 


and Marion О. MekInney, 1951, 
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FigurRE 1.—The stability system of ares. Arrows indicete positive directions of moments, 
forces, and angles, This system of axes is defined as an orthogonal system having the origin 
at the center of gravity and In which the Z-axis is іп the plane of symmetry and perpendicu- 
Jar to the relative wind, the X-axis is in the plane of symmetry and perpendicular to the 
Z-axis, and the Y-arls is perpendicular to the plane of symmetry. Ata constant angle of 
attack, these axes are fixed [n the airplane, 
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average fuselage height at wing root, 
feet | 

average fuselage width at wing root, 
feet 

vertical distance of quarter chord of 
wing root chord from fuselage 
center line, positive downward, feet 

nondimensional time parameter based 
on span (Vt/b) 

longitudinal distance rearward from 
airplane center of gravity to wing 
aerodynamic center, feet 

longitudinal distance from leading 
edge of vertical tail chord to hori- 
zontal tail aerodynamic center, feet 
(see fig. 6) 

vertical distance from horizontal tail 
to base of vertical tail, feet (see 


fig. 6) | 
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height of center of pressure of vertical 
tail above longitudinal stability axis; 
values of 2 must be calculated for 
each angle of attack, feet 

aspect ratio 

sweepback of wing quarter-chord line, 
degrees 

taper ratio (Tip chord/Root chord); 
also, differential operator in Laplace 
transform 

dihedral angle, degrees (sco sketch of 
fig. 9) 

dihedral angle of wing tip, degrees 

time, seconds 

airspeed, feet per second 

radius of gyration about principal 
longitudinal axis of inertia, fect 

radius of gyration about principal 
normal axis of inertia, feet 

radius of gyration about X-axis, feet 
(Jk хо cos? 7-2,2 sin? 1) 

radius of gyration about Z-axis, feet 


(Akzj соз? п Кх, зіп? т 


product-of-inertia factor 
(82,2 — kr.) Sin 9 COS y) 


a - ** -- - 


angle of attack of principal longitudi- 
nal axis of inertia, degrees (sce fig. 2) 

angle of clunb, degrees (sce fig. 2) 

angle of attack of longitudinal body 
axis, degrees (see fig. 2) 

angle between principal longitudinal 
axis of inertia and longitudinal body 
axis, degrees (sce fig. 2) 

air density, slugs per cubic foot 

angle of bank, radians 

angle of yaw, radians 

angle of sideslip, radians 

rolling velocity, radians per second 
(d¢/dt) 
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FIGURE 2—System of axes and angular relationship in fitght. Arrows indicate positive 


P, Pa... P; 
Ar Aa А, а 


direction ot angles. y=a—e. 


yawing velocity, radians per second 
(dy /dt) 
initial angle of bank, radians 
initial angle of yaw, radians 
initial angle of sideslip, radians 
nondimensional initial rolling velocity 
dojde) 
nondimensional initial yawing velocity 
(dy/dc) 
Routh's diseriminant or real part of 
complex root А+ Л 
imaginary part of complex root В 
coefficients of the characteristic bi- 
quadratic equation 
factors of the В, С, and D coefficients 
roots of characteristic biquadratic 
equation 
differential operator (d/de) 
period of the lateral oscillation, seconds 
time to damp to one-half amplitude, 
seconds 
time conversion factor (m/pS V) 
nondimensional time factor (т) 
relative density factor (m/pSb) 
impressed rolling moment, foot-pounds 
impressed yawing moment, foot- 
pounds 
impressed lateral force, pounds 
impressed rolling-moment coefficient 
impressed yawing-moment coefficient 
impressed lateral-force coefficient 
lift coefficient (Lift/¢S) 
drag coefficient (Drag/gS) 
rolling-moment coefficient 
(Rolling moment/qSb) 
yawing-moment coefficient 
(Yawing moment/¢S6) 
lateral-force coefficient 
(Lateral force/¢gS) 
dynamic pressure, pounds per square 


foot ( =P г) 


272483—51——53 
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n 
"чака 
Ca, 
те. 
Cy, 
У. 4 и 
pur Ci, 
е ек, 
аб, _ . 
"=R; 
Cy 
Y= 9 " 
(АС), increment in C4, produced by lift and 
induced-drag forces 
(AC3,), increment in C4, produced by drag 
not associated with lift 
horizontal tail 
а, section lift curve slope 
Subseripts: 
wing wing 
fus fuselage 
tail used to designate vertical tail 
design used to designate design under con- 
sideration 
data used to designate design for which 
force-test data are available 
exp experimental 
V-tail V-tail 
е effective 
H horizontal tai 


CALCULATION OF LATERAL STABILITY AND 


RESPONSE 


Various types of calculations may be performed to indicate 
in some way the stability of an airplane or the response to 
gust disturbances and control manipulations. The calcu- 


lations most commonly made are calculations of time his- 
tories of disturbed motions, period and damping of the free 
motions, and spiral and oscillatory stability boundaries (lines 
of neutral damping of the spiral mode and of the lateral 
oscillations). Step-by-step procedures for performing these 
types of calculations are explained in the text and derivations 
and additional pertinent material are presented in appendixes 
A to D. | 

The period and damping calculations are the easiest of 
the three types to perform. For this reason, and because the 


dynamic lateral stability of airplanes is at present specified 
in the flying-qualities requirements in terms of the period . 
and damping of the lateral oscillation, period and damping 
calculations are probably the most commonly performed. 
Recent dynamic stability work: has indicated, however, 
that the period and damping characteristics of the free’ 
motions of an airplane are not always a sufficient indication 
of whether the dynamic behavior of an airplane following 
various types of disturbances will be considered satisfactory. 
For this reason the calculation of time histories of the 
motions of airplanes is becoming more common despite 


. the fact that these calculations are fairly laborious. The 


increasing use of automatic computing machines has also 
made the calculation of motions more popular. 

For many years, calculations of stability boundaries were 
the type of calculation most commonly performed. In 
recent years, however, stability boundaries have not been 
considered to give an adequate indication of stability. 
Since boundaries are useful in some cases, however, (for 
example, for quick approximation of the effects of changes 
in dihedral and tail area) the methods of calculating the 
spiral and oscillatory stability boundaries are described 
herein. Lines of constant period and damping of the lateral 
oscillation are related to stability boundaries (lines of 
neutral stability). In some cases these lines of constant 
period and damping may prove more useful than boundaries. 
Since no extensive use has been made of lines of constant 
period and damping, however, the methods of calculating 
these lines (presented in references 8 and 9) are not given 
in the present report. 

The equations and methods of calculation presented in 
the present report deal specifically with the inherent motions 
of airplanes for the case of three degrees of freedom (roll, 
yaw, and sideslip) and linear stability derivatives. In order 
to perform similar calculations for cases involving additional 
degrees of freedom, nonlinear derivatives, or autopilots with 
time lag, special equations are required. The methods and 
equations for treating these cases are presented in references 
10 to 18. Additional degrees of freedom for the case of free 
controls are treated in references 16 to 18 and for the ease 
of fuel sloshing are treated in reference 10. The use of 
nonlinear derivatives in stability calculations is covered ш 
reference 11. Methods of treating the effect of autopilots, 
including the effect of time lag in the autopilot, are presented 
in references 12 to 15 and 19. 

For some cases the effects of aerodynamic time lag are 
important. There are two different sourees of such lag: 
(1) the time required for an aerodynamic impulse to travel 
from one component of the airplane to another (for example, 
the time required for a change in sidewash at the wing to 
reach the tail—a phenomenon commonly referred to as lag 


METHODS FOR CALCULATING LATERAL STABILITY AND ESTIMATING STABILITY DERIVATIVES 


of sidewash); end (2) the time required for the growth and 
decay of the aerodynamic loads on the airplane components. 
For both of these cases the time-lag effects usually become 
increasingly important as the period of the lateral oscillation 
decreases. The effects of the first type of time lag сап be 
accounted for in some cases by modification of the stability 
derivatives. For example, the effect of the lag of sidewash 
on the derivative С, is discussed subsequently under the 
section on "Estimation of Lateral Stability Derivatives." 
In many cases, however, both types of time lag will require 
special stability equations. No general treatment of these 
cases has been published but an indication of the method of 
treatment may be obtained from the treatments of autopilot 
lag in references 13 апа 15. 


CALCULATION OF PERIOD AND DAMPING . 


As pointed out in references 1 and 2, the period and damp- 
ing of the various modes of the lateral motion may be calcu- 
lated from the roots of the characteristic equation 


ANM+BM+CMN+DAFE=0 


by the equations 


p= г 
and 
log, 3 
en „0803, | 


where £ represents а real root X or the real part of a complex 
root A=- № and Г represents the imaginary part of a 
complex root. Negative values of Та represent the time 
required to double amplitude for unstable modes of the 
motion. 


The values of the coefficients А, B, C, D, and E may be 
obtained by the method given in steps 1, 2, and 3 of the sec- 
tion on “Calculation of Motions." If the period and time 
to damp are to be calculated for а number of related cases, 
however, the values of the coefficients A, B, C, D, and Е may 
be more conveniently calculated by & tabular procedure such 
as that shown as table I for making boundary calculations. 


Methods of determining the roots of the biquadratic 
characteristic equation are presented in appendix С. 


CALCULATION OF MOTIONS 


Caleulation of the lateral motions of ап airplane involves 
the integration of three simultaneous differential equations 
(see appendix À) to obtain & general solution in terms of the 
mass and aerodynamic parameters of the airplane. The 
general equations, once obtained, can then be used to obtain 
numerically the motions of any airplane in terms of the 


527 


variation with time of the angles of bank, yaw, and sideslip or 
some function of these angles such as rolling or yawing 
velocity. Various methods, such as those given in references 


20 to 22, are of course available for integrating the differential - 


equations. Since the problems met in airplane dynamics are 
fairly complex, however, many of these methods are not 
suitable because of the difficulties of computation that arise. 
The method given in reference 4 (based on the Heaviside 


operational calculus) is satisfactory for calculating the forced | 


motions following application of external forces or momenis 
but, without modifieation, this method cannot be used to 
calculate the motions resulting from initial displacements 
in bank, yaw, or sideslip or from initiel values of rolling or 
yawing angular velocity. A solution based on the Laplace 
transformation is more satisfactory than that based on the 
Heaviside operational calculus because it permits direct 
calculation of the free motions following any initial condition, 
in addition to calculation of the forced motions following 
application of external forces and moments. "The application 


of the Laplace transformation to the calculation of lateral — 


motions is outlined in appendix B. The material presented 
in this appendix is similar to the work presented in references 
5 and 6 except that the mass and aerodynamic stability 
derivatives have been combined as shown in appendix A 
to reduce the number of arithmetical and algebraic processes 
required in numerical solutions. | " 

The process of calculating the motions is presented as a 
series of simple though lengthy arithmetical апа algebraic 


steps so that an understanding of the calculus involved in | 


solving the differential equations is not required. The 


method es shown is suitable for calculating the motions as. 


variations of ф, у, B, р, and r with time for the case of the 
free motions following initial angular displacements (фо, 
Vo, апа Ва) and angular velocities (Dé), апа (Dy), and for the 
ease of the forced motions resulting from constant impressed 
forces and moments (Le, N., and Р.). These are the cases 
for which motions are usually calculated. It is also possible 
to calculate the motions resulting from impressed forces 
and moments which are arbitrary functions of time by the 
methods explained in references 6 and 7. 


MOTIONS RESULTING FROM INITIAL ANGULAR DISPLACEMENTS AND 
ANGULAR VELOCITIES AND FROM CONSTANT IMPRESSED 
FORCES AND MOMENTS 


The six steps involved in obtaining a specific solution for . 


the lateral motions of an airplane are: 
Step 1: Determine values of the following parameters: 
(а) Mass characteristies 


т, Er, Као 7, апа p 
(b) Geometrie characteristics 


S and 6 


+ 
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(c) Flight conditions 
V, Cr, and v 
(d) Aerodynamic stability derivatives 
Сы, Сав Сур Съ Съ» Су С См and Cy, 


The methods of determining Ше values of the aerodynamic stability derivatives are given in subsequent sections of this report. 
In cases where impressed forces and moments are used as disturbances, determine the values of the factors 


C C, Съ, 


that are appropriate to the particular problem. | | | | 
Step 2: From the known factors, evaluate the following parameters which are the stability derivatives in the form in . 


which they are used in the calculation of motions: 


Kxz | К, Куа m X m | = 








ARa Kr "pSV 825% 
3 Съ пау ; Cs, Ui Съ, 

а С, =з Cs, иа, Cr, 

„тр С. me ет Cs, ET Съ, 


Also, when impressed forces and moments are used, evaluate 


Lara С, a ge C, Ve Cr, 


The values of Kz’, Kz?, and Kyz can be determined from the following expressions: 
КК? cos? 7 TÉR sin? q 
Kz-Kz! соз? --Kr sin? 4 


Куз=(К.2— Ке?) sin 7 COS q 


where 
ky, 
Kx P 
k 
Ky 27 


Step 3: Solve for the values of Ше appropriate ones of tbe following coefficients from equations (1) to (4): 
In all cases solve for the values of A, B, C, D, and E: 

А=1— КК, 

B=P,—Ays 


C= Pyet Pat Psy, Ра — Р, 
C C (1) 
D P, +P, = tan y+ P; 


Е=Р, +P, tan y 
where l 


P, = - — n HE mnp tH КУ, 
Ру pn, — lny . 
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P,—lgn,—lng 
Pi=lpnse— isn, 
Р.= Купе і, 
= К; — пв 
Р,=— Pyet Piy Ра. Pi 
The quantities P, to P; are factors of the coefficients B, C, D, and E which are combinations of terms that occur frequently 
in calculations of motions resulting from initial angular displacements and velocities and which are consequently grouped 


tugether for convenience. 
Calculate the values of aq, a1, . . . а; when solving for the angle of bank ¢ or the rolling velocity p: 


Q= hÁ 
а-<-фВ--(ПфуА . 
а= GC —BsP,- (DeC — дув Kul — n,) (Du Kn, — Hien K 


аз--фо ( Р, e tan v+P:)—WoPs e tan y—BsP,-F (Dé) Psyc— Ps— Kal ya trys) CDU — Psy; + 


P,— Ran,ya—l yg) — (n, ya) J-n (aya ti) y. Ps " 
a= [pP PoP t (D$) P — ID PU e tan y+l{ng—ngy-t+ nye) +nllsy,—ls—lLyp)—yePs 
ав-(--ла пай) Е tan y i 
Calculate the values of bo, бъ... bs when solving for angle of yaw y ог the yawing velocity r: 
bo= poA 
| бү:=%В-+ (ОА s | 
i — ————— 
= ФР, (Р, E+ P.) ВРО Pur, + Kal n9 D CP — 
Kingyetl ya) 3- (Kaye n5) —n«(U, tH y) — y. Ps " 
һ={—фР,+%Р,—(РфРу+(Ру»РД E+ ln nry) ndis) P. 
=(lng— ndg) es 
Calculate the values of e, є, . . . С, when solving for the angle of sideslip B: 
со ae 
шеф А С Серда C = tan у+8Р, (Оф Ау, -ОфАф-04уг4 
C4— фе Р, Ery Vo Pi с. tan у+ В. Р. (Оф) Е СК, Гу +-Къ!, Enn, EUG -- пу, |+ 
(Dy), | A ten YT Kin,y,—Kin,y—ly +1, — (E жни у. + Esby)t n —1—Riy,)t y. Fi 
| (4) 


аР, Р, SE tony (D) ( — RE, == Ca tan у-Ел„ С2 tan у+ КА, Cin, E) 


(Dy) ( Kings Се tan Y—i es tan y — Kn, > 2 сън, рия (л, Yr-—Np—Nr Yy tA get tan Y) T. 


Са 


п. (—1, Yrtit y 4 — H3 = fe tan y H- yc P4 


ect (п, Stan y—n,; AN (o5 Cy — —; CS tan Y) 
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Step 4: Solve for the roots X;, Аз, №, and А; of the biquad- 
ratic equation 


АМ - B 4- СМУ РА Е —0 (5) 


where the values of the coefficients А, B, . . . etc. were 
given by the solution of equations (1). Methods of deter- 
mining the roots of the biquadratic equation are given in 
appendix C. 


Step 5: Use the coefficients obtained from equations (1) 
to (4) and the roots of equation (5) to solve for the following 
coefficients: 


Calculate the values of the factors А, 43, . . . 4, when 
solving for the angle of bank ¢ or the rolling velocity p: 


ub = а ал аз ал? GA, Нав . 
ВАА, 5--5ВА, HAO BDA +E, 


P ааа Наз 4- а + 4434 Нав, 
? 6AXST-5BAS--4CXS-E3DX£--2EX, 


iam аз ам Нам ам + Gada t а» 
s BAMI +5 Bist t+ 4 Ons 3-3 DA +2E M 


A= Aore+ 244 + 1 Ag+ 04334. 3- 044, - as 
t 6AAX£54-5BXAH-4C0X38-3DA-2EX, 


(6) 





А,= 32 


Ac (а‹—а, 2) 


Calculate the values of the factors Bj, B5, . . . В, when 
solving for the angle of yaw y or Ше yawing velocity г: 


B= baht В.А DaM tH Бад 2 b+ b; 
GAN +5 BA +4CK + 3DA ФЕ», 


B ENS hs bAt- ,%3-- Б.А? 6 b, : 
7 0AXS--5BX,--4CA3--3DAS--2EX, 
в. Dou ВА bA t дама E bAt Da 
5 6AX,42-5BX,-- ACA) + BDAY -2£EX, 


B bods t 6, ba) Бл + bidet b, , 
+ 6AXPT-5BAT-40AT3DA22EX, 


(7) 


Calculate the values of the factors С, Ci, . . . C, when 
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solving for the angle of sideslip f: 


С tM Ee eut Кам ea 
1 бАХ8-|- SB rttr ЗРАРР РАМ 


Ах БВС + 8D Ae 2 Е 


RE Corgi C4 A3 F 02,3 H- 634? F 644 (8) 
BAA +S Bg EFACAS4-3DX--2EX, 


Дева Co FH eM е 4-с  H- с А, 
BAM HE SBN IONE TDN ZEN 


"QA 
B 

If equation (5) has conjugate complex roots, the values of 
the coefBeients (equations (6) to (8)) corresponding to these 
roots wil be conjugate complex. In order to facilitate 
treatment of this case it is convenient to establish some 
special notation. This special notation is explained in 
appendix D. 

Step 6: The equations of motion are written in different 
form depending upon the roots of equation (5). Н the 
characteristic equation has four real roots М, №, №, and X, 
the general form of the equations of motion is used as follows: 


$— A e™1 + Aht Ае%-|- Ayers Aga + uA, 
y — Bye™1-+ Bye”: -I- Byes Be™+ Boo ВБ, 


В = Се” + Cae + Се" 4- Следа + С 
(9) 
р=> (Ален ме? Ahe Дочев Ар 


К =< (Byd,e"1-+ Borges + Bare + Bye +В) | 


If, as is generally the case, equation (5) has two complex 
roots and two real roots (4-2, R— i, Xs, and X), the equa- 
tions of motion тау be expressed as 


p= Kae"? cos (ol -+ wa) + Age s+ A 074 d- а +A, 
Y= Къе" cos (oI +wp) 4- Bz” + B,e4-- Ва 4- B, 
B — Kece" cos (of ек) + Cses + Ce" С, 

p= | КА +1 26°F cos (of d-o4--tan^! 5)* 


(10) 
А Ме?з + АХ (е ча + А. 


r= | ШТЕР eF соз (oI +us+ tan^! nyt 


Br зе? а + В.А," + В, | 
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where 
K= R+ L^ to 4 = tan^! 2 
А 
Ka—24 Ry! + Is? овес Вап 1 a (10а) 
Ко ВЕТ? wetan 32 
ic 


and Ra and 7, are defined in appendix D. 


If there are four complex roots (А-Л, R—It, R' - I, and 
р--Р1), the equations are 


$ — Ке" cos (cT --o4)--K 4 е cos (ol? tog + 
Ас As 

y= A ge** cos («Г twat Kee cos («Г -ов)-+ 
Во + B, | 

В=К—е°® cos (eI ес) | Ке є cos («Г + ес + C, 

p=- | EEF en cos (oT -+o4+ten~ 5) An (11) 


ААТ? er? cos («T -- «4! Нат”! г) | . 
r= | KEPT е cos (Тор апт 5) + Bet 


=I =e | | 


where 











"RE NES EE DS М 

K= 24 R, 74+,” wa —tan^! a 

А 

___ г 
K,’=2yR;7+ B ов —tan^! Ir (1 Та) 

B 

Ko —24 Re+ IU ос —tan^ ре 

C 


The coefficients Ky, Кв, Ke, wa, ов, and ос are defined in equa- 
tions (102) and Ay, Г, a’, and I,’ are defined in appendix D. 
Solve the appropriate ones of these equations of motion 
(equations (9), (10), or (11)) by substituting values of the 
nondimensional time factor с in the equations and solving 
for 6, у, B, p, or г. 
MOTIONS RESULTING FROM ARBITRARY DISTURBANCES 


The motions resulting from arbitrary forcing functions 
can be obtained from the motions resulting from constant 
impressed forces and moments by the methods explained in 
references 6 and 7. 

A very useful method of obtaining the motion resulting 
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from various abrupt gust and control disturbances is given. 
by Jones in reference 7. In this report it is pointed out that, 
although the component motions of an airplane must be 


calculated simultaneously (that is, by simultaneous differ- - 


ential equations), the effects of component disturbances may 
by the principle of superposition be calculated separately 
and later added in any desired proportion. "Thus, if a given 
rolling moment causes & 20? bank in 1 second and if & given 


yawing moment causes а 5° bank in 1 second, the combined . 


effect of both acting simultaneously will be а 25° bank in 1 
second. Jones also points out & somewhat similar fact with 
regard to the effects of disturbances that are not applied 
simultaneously. This fact is that, if a given disturbance 
which arises at the time & is later augmented, the effect 
of the increment of disturbanee will run its course inde- 
pendently of the effect of the original disturbance. For 
example, in a problem involving the correction for a gust 
disturbance by & manipulation of the control, the motion 
produced by the gust disturbance can be caleuleted inde- 
pendently апа the motion caused by the assumed corrective 
control manipulation can be added to it at any desired point. 
This example is illustrated graphically in figure 3. 

The principle of superposition may be applied analytically 
as well as graphically. The analytical application which 
makes use of Carson’s integral or Duhamel’s integral is 
described in references 7 and 23. This method is useful for 
calculating the motions resulting from impressed forces and 
moments which are arbitrary functions of time. By applica- 
tion of these methods, the solutions for constant impressed 
forces and moments can be used to obtain new solutions for 
any arbitrary variation of impressed forces and moments 
with time which can be expressed by a mathematical formula. 
Some simple variations of impressed forces and moments 
with time and their Laplace transforms are given in reference 


—_— = 


6. The transforms for any other function for which trans- 


forms have been worked out may be found in tables of 
Laplace transforms. 


CALCULATION OF STABILIFY BOUNDARIES 
OSCILLATORY STABILITY BOUNDARIES 


As pointed out in the preceding section of this report, the 
degree of stability of the uncontrolled motions of an airplane 
is indicated by roots of the characteristic equation 


AX*H- BN HT CN D+ E—0 


For stability the real roots or the real part of the complex 
roots of the characteristic equation must be negative. A 
useful discriminant for determining some of the character- 
istics of the roots in stability work is Routh’s discriminant 
В (R=BCD—AP—B'E). The use of this discriminant in 
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dynamic stability analyses has been pointed out in many 
reports, for example, references 1, 2, 3, 5, 21, and 24. Routh 
has shown (reference 20) that, if R and the coefficient Е are 
finite, the necessary апа sufficient conditions that the real 
roots and the real parts of the complex roots should be 
negative are that every coefficient of the biquadratie and 
also В should have the same sign. Routh also showed that 
when Л-<-0 and В and D have Ше same sign there are a pair 
of complex roots with the real parts zero. Since the value of 
the real part of a complex root indicates the stability of an 
oscillatory mode of the motion of an airplane, the lateral 
oscillation is neutrally stable when R=0 and the coefficients 
B and D have the same sign. Oscillatory stability boundaries 
can be determined, therefore, by solving the equation R=0 
and checking to determine whether the signs of В and D are 
the same. 





(a) Gust disturbance. 
(b) Control manipulation. 
(c) Gust disturbance and corrective contró) manipulation. 
FIGURE 3,--ШизизНоп of superposition of motions to determine effect of arbitrary 
disturbances, 
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Since two of the most important stability derivatives af- 
fecting lateral stability are the directional stability derivative 
Crs and the effective dihedral derivative С boundaries for 
neutral oscillatory stability are usually calculated as a func- 
tion of these two derivatives as illustrated in figure 4, These 
calculations are generally carried out by the method shown 
in table I. This table contains a numerical example and 
step-by-step instructions for using the table. The results of 
this numerical example are plotted in figure 4. The pro- 
cedure illustrated in table I is first to assume values of the 
independent variable C,, to cover the range for which the 
boundary is required. The values of all the other mass and 
aerodynamic stability derivatives except C, are then esti- 
mated. The value of C,, is generally assumed to have been 
varied by varying the size of the vertical tail and, con- 
sequently, the tail contribution to each of the other stabilit y 
derivatives varies as Ca, is varied. The values of the coef- 
ficients A, В, C, D, and Е and then 7? are caleulated as 
functions of В: 


m 
“oR Ch, 
The values of /; corresponding to the assumed values of On, 
for the condition of neutral oscillatory stability are пе 
obtained by solving Ше expression /?=0 which is a quadratie 
in в that is of the form 


alg. + ова w= 0 


Finally, the values of С, corresponding to the assumed 
values of.C,, are obtained from the values of dp. 

The values of fs which satisfy the expression Л = 0 must be 
checked to determine whether they satisfy the other соп- 
dition for neutral oscillatory stability—that the sign of the 
coefficients B and D must be the same. This check сар be 
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Бахов 4.—Lateral-stability boundaries calculated in table I. С; g Was the dependent vari- 
able. Cpg was the Independent variable. Сад was actually varied by changing Cy, a 
ta. d 


Varying 


C, Р їп this manner caused changes їп the tall contribution to all the other deriva- 
tives, 
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performed readily by substituting the values of fg which 
satisfy Д--0 into. the expression for D which is a linear 
equation of the form 


Р =з 


Thus, the sign of D is determined. The sign of В is a con- 
stant for any given value of Ca, and is almost invariably 
positive since the three predominant terms of B contain the 
derivatives Су, Са» and Cy, which in all practical cases 
contribute & positive increment to the value of B. 

Since two values of Cj, satisfy the condition Ё=0 for each 
value of Ca, the ?=0 curve has two branches. Аз pointed 
out in reference 24, one of the branches of the R=O curve 
generally represents an oscillatory stability boundary and 
the other braneh represents а line of numerically equal real 
roots with opposite signs. (See fig. 4.) If neither of the 
values of С. which satisfy the expression //—0 for а, par- 
ticular value of С. is found to represent a point of neutral 
oscillatory stability, the lateral motion has no oscillatory 
mode for that value of C., Н both of the values of С, 
which satisfy the expression Ё=0 are found to represent 
points of neutral oscillatory stability, the latera! motion 
has two oscillatory modes. Їп {из case, since the boundary 
D=0 represents the line of infinite period, the branch of the 
И-<0 boundary which lies close to the D—0 boundary is 
usually the boundary for neutral stability of the longer period 
of the two oscillatory modes. 
significance of the stability boundaries and the regions formed 
Бу these boundaries 15 given in reference 24. 

In ealeulating stability boundaries for & specific airplane 
a complete solution such as that explained in the preceding 
paragraphs should be made. For general studies of stability, 


however, &p реше oscillatory stability boundaries may 
be calculated much more simply by the methods shown in 


reference 24. 

Аз pointed out previously, methods of calculating lines of 
constant period and damping of the lateral oscillation are 
presented in references 8 and 9. 


SPIRAL STABILITY BOUNDARIES 


Spiral stability boundaries, like oscillatory stability bound- 
aries, are usually determined as a function of the directional 
stability derivative C,, and the effective dihedral derivative 
Са, as illustrated in figure 4. Аз pointed out in reference 1, 
neutral spiral stability occurs when the Е coefficient of the 
characteristic equation is zero (Е--0). А spiral stability 
boundary can be easily obtained from this relation. If 
expressions for Е (in terms of Їз) corresponding to several 
values of Cx, have already been obtained in the process of 
calculating an oscillatory stability boundary, the equations 
formed by setting these expressions for Е equal to zero can 
be solved for the values of ls (and hence C) corresponding 
to the assumed values of C,,. If the values of Е have not 
already been obtained in the process of calculating an oscil- 
latory stability boundary, а spiral stability boundary for the 
level-flight condition (y=0) can be calculated simply from 


112183--94----84 


A detailed discussion of the- 


the equation 


Values of C,, are assumed within the range for which the 
boundary is required. The values of C; and CX corre- 
sponding to each value of C,, are then determined, The 
tail contributions to these doricauiees generally vary with 
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Cx, Since Cr, is usually assumed to be varied by changing | " І 


the size of the vertical tail. 


ESTIMATION OF LATERAL STABILITY DERIVATIVES 
GENERAL REMARKS 


Methods of estimating the lateral stability derivatives have 
been presented in numerous publications but no single report 
has contained information for estimating the contribution of 
all principal airplane components to all the derivatives for 
airplanes having any sweep angle or aspect ratio. In the 


~ + ee 


present report, ап approach to such a presentation is made 


by the coordination of and reference to existing estimation 
methods, by reference to publications containing data which 
should be useful in making estimates, and by the suggestion 
In some cases of simple new empirical formulas. Detailed 
estimation methods аге presented for low-subsonic-speed 
conditions but only а brief discussion and а list of references 
аге given for transonic- and supersonic-speed conditions. 
In general, the estimation methods presented should be ex- 
pected to yield only fairly accurate values suitable for making 
first approximations of dynamic stability. "This limitation 
applies especially to the cases in which the derivatives are 
based completely on theory. 


For convenience, the references that should be useful in. 


estimating the stability derivatives are presented in table II. 
The references are grouped according to the speed range 


covered (subsonic or supersonic) and according to the deriva-. 


tives presented in each report. The references for the sub- 
sonic case (references 1 and 25 to 97) are further divided into 


two groups—one including pd aides which contain estimation 
methods and the other including reports which contain 


experimental data that should be useful in making estimates 
of derivatives. The references for the supersonic case (геЃег- 


ences 98 to 118) are subdivided according to wing plan form. _ 


The following sections covering the estimation of the nine 


stability derivatives are divided into three groups according 


-.---- . —— 


to the type of derivative—sideslip derivatives (Су Cag, Cig), 
rolling derivatives (Cx, Cip Crp), and yawing derivatives 
(С. Съ, Съ,). The derivatives Съ, and Cy, have usually 
been neglected in making dynamic lateral stability caleula- 
tions because theory indicated that for unswept wings Cy, 
and Cy, were zero. Recent experimental data, however, 
have indicated that both swept and unswept wings produce 


measurable values of these derivatives (references 25, 60, 


and 89). Since the vertical tail contributes to Cy, and Cy. 
it appears desirable to estimate these derivatives ino to use 


* 
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them in the calculations of stability unless it is established 
that for the case in question the effects of Су, and Су, on 
stability are negligible. For these two derivatives, only the 
effect of the wing and vertical tail need to be considered. 
The methods of estimating the rolling and yawing deriva- 
tives presented herein were obtained from theoretical treat- 
ments based on the assumption of steady rolling and yawing 
and from experimental data obtained principally from tests 
made under conditions of steady rolling and yawing. The 
only information that applies directly to the oscillatory 
case is а limited amount of data on C, obtained by oscillation 
techniques. When calculations ате made in which the 
oscillatory mode is the subject of interest, some consideration 
should be given to correcting the derivatives based on steady 
rolling ог yawing to account for differences in the derivatives 
that are likely to exist as & result of differences between 
the oscillatory motion and the steady rolling and yawing 
motion. For example, the data of reference 85 have indi- 
cated that, for flap-extended or power-on conditions, fairly 
large differences might exist between the values’ of the tail 
contribution to C», for the steady yawing and yawing oscil- 
lation cases. At present little information is available for 
correcting the values of С, for the steady yawing case to 


apply to the oscillatory case and, unfortunately, little or 
no information is available for correcting the other stability 
derivatives. 


Since most wind-tunnel forco-test data that are likely to 
be used in making estimates of the stability derivatives 
are probably for much lower Reynolds numbers than those 
for the full-scale airplane, some adjustments to the data 
are usually required to account for the differences in Reynolds 
number. The effects of Reynolds number should be con- 
sidered in the cases of all the derivatives, especially those 
which are estimated by methods that involve the use of 
force-test, data. Methods of correcting for Reynolds num- 
ber eftects for some of the derivatives are discussed in the 
following sections which cover the estimation procedures. 
In the cases where the Reynolds number effects are not 
discussed, 1, can be assumed that any abrupt variation in 
the derivatives near the stall for low-scale data will also 
be present for the full-scale airplane but will probably occur 
аі а higher lift coefficient because of the higher maximum 
lift coefficient. of the airplane. Án indication of the lift- 
coefficient range over which the theory may not be expected 
to give reliable values of stability derivatives for the full- 
scale airplane can be obtained from large-scale drag data. 
The analysis of reference 89 indicates that the variation 
of the derivatives with lift coefficient. is different, from the 
theoretical variation at lift coefficients above that at which 
the drag due to lift increases abruptly from the ideal value 
C, T. А. 

The effects of Mach number and power are not treated 
in the sections on the individual derivatives but ere dis- 
cussed briefly in separate sections. A detailed treatment 
of these effects, including design formulas апа charts, was 
considered beyond the scope of this report. 
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THE SIDESLIP DERIVATIVES Сур Са, С, 


No satisfactory purely theoretical methods have yet been 
developed for obtaining accurate estimates of the sideslip 
derivatives Су, Cap and Ci, for a complete airplane, pri- 
manly because of large interference effects between the 
various airplane components and because of large, and often 
unpredictable, variations of the derivatives with angle of 
attack. Fortunately, these derivatives can be obtained 
from conventional wind-tunnel force-test data. Such experi- 
mental data are essential to the accurate determination of 
sideslip derivatives. It is, of course, highly desirable to 
have force-test data for the exact airplane design under 
consideration, but reasonably accurate estimates can usually 
be made by correcting the force-test data for а generally 
similar design. The methods of correcting the force-test 
data on а similar design for use in the case under consideration 
are covered in the followmg sections. In the formulas 
presented, the subscript word “design” is used to designate 
the design under consideration and the subscript word 
"data" is used to designate the similar design for which . 
force-test data are available. 

Force-test data should be used to determine the effect on 
the sideslip derivatives of such airplane components as 
leading-edge high-lift devices, stall-control devices, trailing- 
edge flaps, nacelles, external stores. canopies, and dorsal and 
ventral fins. The effect of leading-edge high-lift devices is 
usually merely to extend to a higher hift coefficient the same 
variation of the derivative with lift coefficient as for the plain 
wing. Trailing-edge flaps often have large effects on the 
contributions of both the wing and the vertical tail to tho 
sideslip derivatives (references 40 and 71); and since these 
effects are not easily estimated, it appears that in these cases 
use of force-test data is essential. The addition of nacelles 
and externa! stores generally has been found to decrease the 
directional stability factor Ca, slightly. The results of а 
limited amount of research to determine the effect on the 
sideslip derivatives of the size and shape of canopies has 
been reported in references 49 and 75 but these results are 
inadequate for making accurate predictions of the effects 
of canopies. The effects on the sideslip derivatives of dorsal 
and ventral fins are usually small at the small and moderate 
angles of yaw that are generally considered in stability 
calculations. (See references 48 and 73.) 

Cy, 

In estimates of the lateral force due to sideslip derivative 
Су force-test data for the design under consideration 
should be used whenever possible. If such data are not 
available, data for a similar design can he used and corrected 
as follows: 


Wing-fuselage.—Since the wing-fuselage contribution to 
Cy, is usually relatively small compared with that of the 
vertical tail, great accuracy is not required in estimating this 
factor. This contribution may be estimated as follows: 

(1) Wing: If the wings of the two designs are generally 
similar, the difference in Cs, iss can be considered. negligible 
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and no correction is necessary. The theory of reference 25 
does not appear to be suiteble for use in estimating Се, 


(2) Fuselage: If the two fuselages аге similar in Ерна 
Ше difference in Cry can probably be estimated satisfac- 


torily by correcting for the difference in the relative size of 
the fuselage and wing for the two airplanes. It appears, 
however, from table X of reference 71 unlikely that a reliable 
prediction of Cys, can be made directly from the geometry 
of the fuselage. "Some additional data on Cy. are pre- 


sented in reference 79. Experimental data from other 
investigations have shown that differences in fuselage cross 
section can cause very large differences in the variation of 
Cr. with angle of attack. For example, in the case of & 


flat. fuselage with the major cross-sectional axis horizontal, 
the sign of Cy, has been found to reverse at moderate and 


high angles of attack. Force-test data are essential for 
making estimates m such cases. 

(3) Wing-fuselage interference: For low-wing or high-wing 
configurations, wing-fuselage interference causes the value 
of Cy, to be greater than that obtained by adding the con- 
tributions of the wing and fuselage. (See references 37 and 
40.) If the vertical location of the wing on the fuselage is 
generally similar for the two designs, however, any correction 
for a difference in this interference factor can be neglected. 

Vertical tail Accurate estimates of Cr st are necessary 


beeause this factor is used to estimate the tail contribution 
to several other derivatives. This factor is especially im- 
portant at low angles of attack because in this case the tail 
contribution is often much greater than the wing-fuselage 
contribution to all derivatives except С,. For this reason 
it is highly desirable to have tail-off and tail-on force-test 
data for the design under consideration or for a very similar 
design. Corrections to the data for а similar design can be 
made as follows: 

(1) Correction for differences in wing area, tail area, and 
tail lift-curve slope can be made by the following formula: 


(C ы Ө 
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The value of C, сап be obtained from figures 5 and 6 


which аге based on the theory of reference 34 апа on the 
theory and data of references 28 and 35, respectively. The 
ehart of figure 6 can be used to estimate the change in the 
effective aspect ratio of the vertical tail caused by the end- 
plate effect of the horizontal tail. It should be emphasized 
that for the best accuracy the charts in figures 5 and 6 should 
be used in conjunction with formula (13) for correcting 
existing force-test data and not for making a direct estimate 
of Сү, 


(13) 


(2) In the ease of V-tails, the ГЕТЕА for Cons сап be 
made 85 follows: 
(КС, Sv. Hats sin" Т) ае Sante. 


( ард. (С ZONES DG (KC, S v-tati sin? Г) ыы 


Sura desiga 


(14) 


where Ше terms С. , Г, and К are the same as given in 
reference 30 and are defined as follows: 


Съ 

P plane normal to the chord plane of each tail panel 

г dihedral angle of tail surface measured from .X Y-plane 
of the tail to each tail panel, degrees 

K ratio of sum of lifts obtained by equal and opposite 
changes in angle of attack of two semispans of tail 
to lifts obtained by an equal change in angle of 
attack for the complete tail 


Values of the term K, which are usually about 0.7, can be 
obtained from reference 30. 
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slope of the tail lift curve in pitch measured in the — 


(3) Since large differences in sidewash апа dynamic pres- | | 


sure &t the tail can be caused by differences in wing plan 
form and wing location, use of experimental data for the 
specific design or at least for а design which has а closely 
similar wing-fuselage combination and vertical tail location 
1 extremely desirable. No methods are available which 


—— a——mÁ— 


permit accurate predictions of sidewash at the tail, but the __ 


experimental data of references 40, 50, and 71 can be used to 
obtein some indication of the variation їп sidewash with 
vertical location of an unswept wing on а fuselage and Ше © 
experimental data of references 36 and 79 provide additional 
information on sidewash at the tail Other experimental 
data indicate that the sidewash fields produced by highly 
swept, low-aspect-ratio wings or by fuselages of flat cross 
section can sometimes be strong enough at high angles of 


attack to reverse the effectiveness of а conventionally located — 


vertical tail surface. Until а reliable method is developed 
for predicting these large sidewash effects, force-test data 
appear to be the only means by which satisfactory estimates 
of Cs, и сап be obtained. 


Although attempts h&ve been made to develop methods 
for estimating the yawing moment due to sideslip (static 
directional stability) derivative Ca, (for example, references 
70 and 71), no-reliable method has yet been obtained. The 
use of force-test data therefore seems imperative. 

Force-test data for the design under consideration should 
be used if available. If such data are not available, use data 
for a similar design &nd correct as explained in the sections to 
follow. 


Wing-fuselage.—The corrections for the wing-fuselage 
contributions are: 


(1) Correction for wing: From figure 7 (taken from refer- 


ence 25) the values of (C 


(С. а for the design under 
consideration and for the 4 


esign for which test data are 


available can be determined. The effect of differences in | 


taper ratio can be neglected. (See references 61 and 68.) 
The difference between these values of C,,/C,? should then 
be added (with proper regard for sign) to tlie experimental 
data for the complete model. 

(2) Correction for fuselage: The formula 


C. ——13 e Fuselage volume (15) 
Brus w 
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can be used to calculate Ше Ca, of the fuselage (per radian) · Formula (15) does not include tho effect of fineness ratio 
for the design under consideration and for the similar design | and should not be used for fineness ratios less than 4. This 
for which force-test data are available. The differences formula is an approximate empirical expression which 


between these two values can then be added (with proper | | i. 
regard for е to the force-test data for Ше complete model. впоша покра оаа “to estimgtg: tie зао Mid directly 
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FIGURE §.—Variation of lift-curve slope with aspect ratio, taper ratio, and sweepback for the case of subsonic incompressible flow, ae=mQ.11. "Values from reference 34, 
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FIGURE 8.—-Е бес of horizontal-tafl location on the effective aspect ratio of the vertical taii 
А сад for the case of sabsonfe incompressible low. a=0° Taken from reference 35. 
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but should only be used as indicated to determine а correc- 
tion for force-test date. This correction method should 
not be used in the eases of high angles of attack when there 


are large differences in fuselage configuration: Force-test 


data are essential in such cases. 


(3) Correction for vertical location of the wing: If the B 
designs are generally similar, the correction for the vertical | 


location of the wing on the fuselage сап be neglected. (See 
references 37 and 40.) 


(4) Correction for center-of-gravity position: If the center- 


of-gravity position for the design under consideration is 


appreciably different from that for the design for which — 
force-test data are available, the value of Ca, for the wing- | 


fuselage combination can be corrected by multiplying the 
value of Cy, for the wing-fuselage combination by the 


distance between center-of-gravity positions (expressed in 


wing spans), 
Vertical tail—Corrections to Cap, 


Съ, апа tail length 2/6 can be Hade by the following 


f lail 
formula: 


| 
Or, b 
(Cs, НИ m (C...) date Crne Jusa (16) 


( Ува b data 


The contribution of wing-tip fins to Cy, is treated in 
references 67, 72, and 87. 
1, | 
In estimates of the rolling moment due to sideslip (effective 
dihedral) Ci, force-test data for the design under considera- 
tion should be used. If such data are not available, data 


for a similar design can be used and corrected by the methods 


that follow. 
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Fic Bg 7.—Varlatlon of С. [Cr with aspect ratio and sweep for the case of subsonie incompressible Йот. 1.0; iT. Taken from reference 25, 
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Wing-fuselage.— The corrections for wing-fuselage con- 
tributions aro: 

(1) Correction for wing: From figure 8 (based on reference 
25) the theoretical values of C,,/C, for the design under con- 
sideration and for the design for which data are available can 
be determined. The difference between these two theoretical 
values can then be added (with proper regard for sign) to the 
experimental data. Consideration should be given to scale 
effect, airfoil section, and surface roughness on the value of 
Са for highly swept wings. The lift coefficient at which the 


ОИ variation of С, with lift coefficient departs 


from theory is greatest at high Reynolds numbers and for 
smooth wings with round leading edges. For wings with 
rough surfaces or sharp leading edges the effects of Reynolds 
number on C; are usually small and low-scale wind-tunnel 


data can be used. For airplanes having very smooth swept- 
back wings with rounded leading edges, however, some cor- 
rection should be made for scale cffect when estimations are 
made from low-scale wind-tunnel data. Since no rational 
method has been developed for making such corrections, it is 


suggested that, for lift coefficients higher than that at which | 


the experimenta! data departs from the theory, an average 
of the theoretical &nd low-scale experimental values be used. 
Conservative dynamic stability results will usually be ob- 
tained if the uncorrected theoretical values of Ci, аго used 
because these values are ordinarily greater (more negative) 
than measured values and because the larger negative values 
of C, usually tend to decrease the dynamic lateral stability. 

(2) Correction for wing dihedral: The effect of dihedral 
on Си is treated in.references 29, 40, 52, 59, 68, and 81. 
Correction for the difference in dihedral between the two 
designs can be made by multiplying the incremental geo- 
metric dihedral angle (in degrees) by the factor Cy, obtained 


from figure 9. А plot of Cis, against aspect ratio for taper 


ratios of 1.0, 0.5, and 0.25 (obtained from references 59 and 
68) and а formula from reference 59 for correcting for sweep 
are presented in the upper portion of figure 9. The lower 
chart and formula in figure 9 (developed from reference 68) 
should be used in addition to the upper chart and formula 
of figure 9 to estimate the values of Си for the case of a wing 
with partial-span dihedral. Although this chart and formula 
apply directly only to wings with one dihedral break, they 
can be used to estimate the Cry, for wings with two or more 
dihedral breaks by the method ‘described in reference 68. 
The effect of drooped wing tips and of wing-tip end-plates 
on Съда, should be determined by experimental data since 
no reliable estimation procedure for these effects is available. 
(See reference 67.) 

(3) Correction for wing-fusclage interference: Although the 
contribution of Ше fuselage alone to С, is usually negligible, 
the interference between the wing and fuselage can greatly 
alter the value of C;, of the wing. This interference is such 
that а high location of the wing on the fuselage gives more 
positive effective dihedral (higher —C,) and а low wing 
location gives less positive dihedral than a midwing position. 


This effect is treated theoretically in reference 69 and has 
been studied experimentally in references 37 and 39 tv 43. 
The following simplified expression for estimating the inere- 
ment in C, caused by wing-fuselage interference has been 
developed from the relationships presented in reference 69 
and in other sources: 


AC 1.254 т "OBS U m (17) 


This expression has been found to give reasonably good 
agreement with experimental data for a variety of config- 
urations. It is suggested that values of АС, be calculated 


from this equation for both the design under consider- 
ation and for the design for which force-test data are avail- 
able. The difference between these values can then be 
added (with the proper regard for sign) to the force-test data. 


Vertical tail.—The value of Cryo determined from forec- 


test data on a similar design can "be corrected as follows io 
obtain eos for the design under consideration: 


(€. 
| идда ONE ON (18) 


That b data. 


The results of reference 35 indicate that Су, can also be 


affected by the location of the horizontal tail V respect. to 
the vertical tail. If the two designs have approximately tho 
same horizontal tail size and location, however, this effect 
can be neglected. 


The value of C 
following empirical formula: 


ч for a V-tail can be estimated from the 


(Cu, ei ш — ( C. pasi Jas 


Съ, 
чан (бул 420и sin T) 
b в Г : 








"uH бү ый (Oeics 420. sin r) 
6 sin Г data 


where 6y...4 is the developed (not projected) span of tho 
V-tail, гуми 18 the vertical distance from the center of 
gravity to the chord of the V-tail (positive up), and Г is 
the dibedral angle of the Y-tail. Мого information on 
V-tails can be found in references 30, 62, and 63. 

In the case of а vertical tail located on the wing, there is, 
in addition to the incremental Си produced by Ше tail 
lateral force, an incremental C;, produced by the interference 
effect of the vertical tail on the wing. Since this inter- 
ference effect varies greatly with spanwise and vertical 
position of the tail, it should be determined from force 
tests. Usually the interference is such that & vertical tail 
above the wing gives а negative increment of Ci (positivo 
effective dihedral) and one below the wing gives & positive 
increment of Си. In general, the largest interference effecta 
are obtained with vertical tails &t or near the wing tipa. 
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FIGURE $.—Variation of Сү Р / Cr, with aspect ratio, taper ratio, and sweep for the case of subsonic incompressible ow. Based on method of reference 25. 
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THE ROLLING DERIVATIVES Cn,, Ci, CY, 
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The wing and vertical tail are the only airplane components 
that contribute appreciably to the yawing moment due to 
rolling derivative С... The contributions of the fuselage and 
horizontal tail can usually be neglected. 

Wing.—The contribution of the wing to Ca, can be esti- 
mated from the formula and charts of figure 10 which were 
tuken from reference 89. Although these charts apply 
strictly only to wings having & taper ratio of 1.0, experimental 
data have indicated that they will also provide fairly good 
estimates for taper ratios of 0.50, 0.25, and 0. In the esti- 
mation formula 


_ (Аб), а LOCOS) 


C, = 0, 20 
Ee, (б.у у. (20) 


(Cn). | 


the value of (Со) should be determined, if possible, from 
force-test data obtained at high Reynolds number on the 
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Figure 9,—Ef/fect of dihedral angle on С; 8 for the case of snbsonic incompressible flow (C в 
In radians; Г ір degrees). Taken from references 59 and 68. 


An) ciuem 3 
Ст Ci ро 
Ё ^r (Pis, алые 
(А-4) cos А 
Cis T АА сова (Cip), cn 


where 


REPORT 1098—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


wing under consideration, since low Reynolds number data 
might indicate values of (Cp,), that are too large. For the 
case of smooth wings with а large leading-edge radius and 
low or moderate sweep, it is suggested that ( C»). for the air- 
plane be assumed to be zero at all lift coefficients up to the 
stall. This assumption will result in larger negative values 
of С. than would be estimated from low Reynolds number 
data on (Cp,), and consequently should lead to conservative 
dynamic stability results since an increase in C, in the nega- 
tive direction has been found to cause & veductioni in dy nemic 
stability. The value of (Cp,), for highly swept wings is often 
very large at high lift coefficients, especially for wings with 
rough surfaces, sharp leading edges, or triangular plan form. 
For these cases, values of (Со), determined even from low 
Reynolds number data might lead to reasonably good esti- 
mates of C In all these cases, however, high-seale drag 
data should be used whenever it is available. 

Effect of high-lift devices.—The principal effect of leading- 
edge high-lift devices is to extend to a higher lift coefficient 
the linear variation of C,, with lift coefficient. The formula 
and charts of figure 10 are i direc (ly applicable іо this case. The 


. effect of trailing-edge high-lift devices is not во straight- 


forward, but experimental date have indicated that the 
formula and charts of figure 10 also give reasonably good 
estimates in this case. 





Aspect ratio, 4 
FIGURE 10,— Variation of (AC, )/Cz and (4C, „)„/(Сро)а with aspect ratio for the 
case of subsonic incompressible flow. =1.0. Taken from reference 89, 
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Vertical tail.—' The contribution of an isolated vertical- 
tail surface to C, can be estimated by the following approxi 


mate formula which has also been commonly used to esti- 


mate C uou of a complete airplane: 
; ix 
С LN NS b Сура Q1) 


The values of Сува should be determined from force-test 


data as previously discussed. Instead of the geometric tail 
length //b, it will usually be better to use the effective tail 
length — Cue af? ув 25 determined by force-test data. 


Formula (21) then becomes 


Саа 2 (5) Саа 


Tn the case of Ше conventionally located vertical tail surface, 
however, the rolling wing produces а sidewash at the teil 
which greatly alters the teil contribution to C,,. This side- 
wash causes the values of Capua to be much more negative 


than is indicated by formula (21). This effect is discussed 
more fully in reference 36 in which 15 also presented а method 
for estimating the sidewash. Some preliminary theoretical 
studies have indicated that the effect of the sidewash on 
Cy - varies considerably with tail size and tail location 


and to some extent with wing plan form. А comprehensive 
experimental verification of this theory is planned but as 
yet only а few scattered checks have been obtained. For 
the case of the conventionally located vertical tail surface, 
the following formula has been found to give estimates of 


(218) 


Capui that are in fairly good agreement with experimental 
data: 
[Tez fa 
ст). ] C7 d 
or 
2 үг 
=2 EOM С. (222) 
This formula is based on the assumption that C,, „а 18 Zero 


at 0? angle of attack and varies with angle of attack і in the 
same manner as indicated by formula (21). Formule (22) 
or the method of reference 36 can be used satisfactorily for 
first approximations of Ca, a for most configurations with 


conventionally located еа tails. Гог more accurate 
estimates, especially for configurations having an unusual 
tail size or tail location, experimental data should be used. 

For wings of triangular plan form with vertical tails either 
directly аБоте ог &bove and slightly behind the wing, exper- 
imental data have indicated that neither formula (21) nor 
formula (22) gives an accurate estimate of Cou but that 


an &verage of the values obtained by the two formulas pro- 
vides а fairly good estimate. 


It is obvious that these methods of estimating C, аге 
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only approximate and are open to question in many cases.” 


Experimental and theoretical studies are currently being 
made to provide better methods of estimating См. These 


studies indicate that the sidewash from the fuselage as well 
as that from the wing should be taken into account ш esti- 
mating Cs, 


ДИК а e matis presented herein should be discarded. 


At the present time, however, formula (22) and reference __ 


36 will usually provide much more accurate estimates of 
Cayo than formula (21) which has been in common use up 
until this time. 
МА 
Wing-fuselage.—M ost of the rolling moment due to rolling 
(damping-in-roll derivative) C; of an airplane is produced 
by the wing. The effect of the fuselage can be neglected 
unless the ratio of the diameter of the fuselage to the wing 
span is relatively large (greater than about 0.3). For large 
values of this ratio, the value of C; will be smaller than that 
for the wing alone by an amount that can be estimated from 
a consideration of the area and lateral center of pressure of . 
the wing area included within the fuselage. (See references. 
106, 111, and 115.) 
Wing —The damping in roll of wings has been the subject — 
of many experimental and theoretical investigations. (See __ 


When these methods become available, the . 


references on С, in table ТЕ.) As a result, some methods of ~ 


estimating Ci, have been developed w hich have been found 


to give reasonably good agreement with experimental results. 


The method presented i in reference 81 appears to give suffi- 
ciently accurate estimates of Ci, for zero lift. This method 
is extended in reference 92 to permit the estimation of C, 
over the normal flight range of lift coefficient. Estimation 
charts and formulas from reference 92 are presented in figure 
11. 

High-lift devices —Experimental data have indicated that 


the damping in roll of wings at low and moderate lift соев- _ 


cients 15 not greatly affected by the addition of high-lift 


devices such as trailing-edge flaps, leading-edge. flaps, slats, - 


and slots. 
the lift coefficient at which the sharp decrease in Ci, occurs. 
The charts and formulas of figure 11 can be used to estimate 


the Сп of wings with either full-span or partial-span high-lift 
devices with fair accuracy despite the fact that the method 
is not strictly applicable to partial-span high-lift devices. 
(See reference 92.) 

Wing-tip fuel tanks.—The use of wing-tip fuel tanks 
usually increases the damping in roll of the wing. The 


experimental data of reference 94 for unswept wings indicate __. 


that the magnitude of the increase varies with angle of 
attack and depends upon the wing taper ratio and on the size 
and location of the tanks. Unpublished experimental data 
indicate similar effects of wing-tip tanks on sweptback wings. 


The principal effect of such devices is to increase - 


. The following approximate formula for estimating the incre- 


ment in C; produced by wing-tip tanks at low lift coefficients . 


is based on the limited emount of available experimental deta 
and should not be expected to yield very close quantitative 
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estimates: 


Maximum tank diameter (Kr) 


(AG Wing span 


р Jana = ( б, ) tanks off 


(23) 
where, for symmetrically mounted tip tanks, 


for tanks mounted below the wing tip or forward on the wing 
tip, 
Kr=3 


and for pylon-mounted tip tanks, 
Кт-- 1 


Experimental data for both unswept and swept wings 
indicate that (ACi ааъ usually becomes smaller with in- 
creasing angle of attack and, in some cases, actually reverses 
sign at high angles of attack so that the tanks are decreasing 
rather than increasing the damping in roll. The data of 
reference 94 can be used to obtain an approximate estimate 
of the effect of angle of attack for unswept wings. 

Tail surfaces.—The contribution to C; of conventional 
type horizontal and vertical tail surfaces is usually very small 
and, in most cases, negligible. When ап airplane rolls, the 
wing produces а rotation of flow at the tail surfaces which 
reduces the already small damping moments of the isolated 
surfaces, except in the case of the vertical tail at high angles 
of attack where the tail center of pressure is below the center 
of gravity. 

The contribution of an extremely large horizontal tail to 
C, might not be negligible and can be estimated by multi- 

? 
plying the value of С. for the particular tail plan form 
obtained from the charts and formulas of figure 11 by the 


2 
factor 0.5 x (3) in which the factor 0.5 is included to 


account for the rotation of flow produced by the wing. 
The contribution of an isolated vertical tail surface to С. 
is given by the following approximate formula: 


г 2 
Ci, T (5) CY t (24) 
AsinthecaseofC, из formula can be modified to provide 


ап approximate correction for the effect of the wing on the 
damping in roll of conventionally located vertical tail 


surfaces: 
zN[2 {г 
Ci, a (5) ROM Cr ut (28) 


An analysis of this expression indicates that the value of 
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FIGURE 11.—Charts and formulas for catimating Су. for the case of subsonie incompressible 
flow. Taken from reference 62, 
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Са is negligible at low and moderate angles of attack 
where z/b is positive but that it might be fairly important at 
very high angles of attack where 2/6 is а large negative value. 
Аз in the case of C, , experimental data indicate that, for a 
vertical tail located either directly above or above and 
slightly behind a wing of triangular plan form, the value of 
Ch, ы COD be estimated with better accuracy by an average 
of formulas (24) and (25) than by formula (25) alone. For 


-conventional tail arrangements, however, formula (25) gives 


better correlation with experimental data. 
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Cy: 
Fp 


Wing.—The following formula for the derivative Су, 
(lateral force due to rolling) from reference 89 is based on 
experimental data and is the same as that presented in 
reference 25 except for an additional correction to account 
for tip suction: 


Tte A+cosA | 


~ A+4 cos at (26) 


an AT 
The data of reference 89 show that this formula applies only 
for lift coefficients below that at which the drag factor 


At higher lift coefficients the 


experimental data indicate smaller values of Су, than given 
by formula (26). For these cases an approximation of the 
value of Су can be obtained from the experimental data of 
reference 89. Asin the case of C, , the break in the variation 
of Cy, with lift coefficient should be expected to occur at 
lower lift coefficients for wings having sharp leading edges or 
rough surfaces and for wings tested at low Reynolds numbers. 

Vertical tail—The discussion concerning С, and 


Ch is also applicable to Су, The value PE 


CH { я 
бот begins to increase. 


for an isolated tail surface is given by the formula = 
2 - 
Cr, tail =2 (5) Ст a (2 r) 


This formule can be modified as follows to account approxi- 
mately for the effects of wing sidewash in the case of a con- 
ventionally located vertical tail: 


С" == LE 
"Es il 7 (5 (5) Де Каш 


An average of formulas (27) and (28) can be used for tails 
located either directly &bove or &bove and slightly behind the 


wing. 


(28) 


THE YAWING DERIVATIVES Cn, Cr, Cr, 


' C 
n 
г, 


Wing-fuselage.—1n the past, the contribution of the wing- 
fuselage combination to yawing moment due to yawing 
(damping in yaw) derivative С, has usually been found to 
be small compared with the contribution of the vertical teil. 
The fuselage contribution to the damping in yaw depends, 
of course, on the relative size of the fuselage and wing. In 
the past, the relative size of these components has generally 
been such that the fuselage contribution could be neglected. 
(See references 85 and 86.) For some recent designs which 
have a large fuselage relative to the wing, however, the 


818. 


fuselage contribution to C, is important. In the case of 
fuselages having flat sides or having a flattened cross section - 
with the major axis vertical the fuselage contribution may - 
also be important and some fuselage contribution to Ca, 
should be assumed, especially at high angles of attack. On 


the other hand, experimental data have shown that a | 


flattened cross-section fuselage with the major axis hori- 
zontal can have negative damping in yaw at moderate and 
high angles of attack. - 

The contribution of the wing to C,, can be estimated from 
the formula and charts of figure 12 which were taken from 
reference 25. Values of Cp, for the wing should be estimated _ 
from force-test data. For values of zíc greatly different from “| 
zero, the charts of reference 25 сап be used. The formula and _ 
charts of figure 12 are not considered reliable at'bigh angles 
of attack, especially forswept wings. The use of experimental 
data from the references on (©, listed in table IT is recom- 
mended in this case.  . 

The effect of partial-span inboard flaps on С, can usually 
be neglected. (See reference 85.) The effect of full-span 
trailing-edge ог leading-edge high-lift devices can be esti- . 
mated satisfactorily from the formula and charts of figure 12. . 
Values of Cp, in this ease are, of course, for the wing with | 
the high-lift device installed. 

Vertical tail. — The contribution of а conventional-type 
vertical tail to Са at low and moderate angles of attack can 
be estimated from the formule 


l 2 
C. 2 (F) Сем. Q9) _ 
or, with the effective tail length — Cui al Су», substituted 
for the geometric tail length i/b, _ 
С. cy 
C, == C C вы) (29а) 


The alternative method of estimating C, presented in 
reference 83 wil probably provide better estimates than 
formula (29) in the higher angle-of-attack range. The experi- 
mental values for Са, " presented in reference 85 for power- 
on or flap-down configurations are 30 to 40 percent greater 
than the values predicted by formulas (29) or (29a). These | 
differences &re &ttributed to lag of sidewash effects in the 
free-oscillation tests used in measuring C,,. In estimations of 
Cana for stability calculations, similar Jag of sidewash 
effects should be assumed if the oscillatory mode is of primary 


. importance but no lag of sidewash should be assumed if the 


aperiodic mode is most important. 
Alethods for estimating the С. for wing-tip vertical 
tails are presented in references 72 &nd 85. 
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FIGURE 12,—Charts and formula for estimating Са, for the case of subsonic incompressible Вот, Taken from reference 25, 
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The wing and vertical tail are the only airplane components 
that contribute appreciably to rolling-moment-due-to-yswing 
derivative C; of an airplane. The contributions of the 
fuselage and horizontal tail can usually be neglected. A 
semiempirical method for estimating С, is presented in 
reference 88. This method involves the use of experimental 
data on the parameter C;, to correct the theoretical values of 


DUE given in reference 25 and to estimate the value of 


Iria р 

Wing.—The formula of reference 88 and the charts of 
©, [Cz from reference 25 for estimating Cu, tu &re given in 
figure 13. The values of Си /Сг to be used in the charts can 
be obtained from figure 8. For taper ratios less than 0.25, 
values of C;,/C; and Cul C for & taper ratio of 0.25 can be 
used. The value of C used in the formula should be the 


same as the value of Cr estimated from experimental 
Poing хр 


data by the method indicated in the section on Су. In the 
ease of C,, however, (unlike the case of Сы) conservative 
dynamic stability results will usually be obtained if the 
smaller values of the derivative (based on low-seale experi- 
mentel data) are used instead of the larger (theoretical) 


values. This difference is a result of the fact that either an 
increase in the normally negative value af C;, or a decrease 
in the normally positive value of Су can cause reduction in 
dynamic stability. 
mation procedure shown in figure 13 appears to account 
satisfactorily for the effects of high-lift е wing dihedral, 


cedure is directly applicable to midwing configurations but 
should not be used for high-wi ing or low-wing configurations - 
because changes in wing position produce much greater 

changes in Cj, than in Cj. Work is now being done to 
develop ап estimation method for the effect of wing position 
on Ci, bui until this method is published the following pro- 
cedure is recommended for estimating О; of high-wing and 
low-wing configurations: Adjust the value of Си to cor- . 
respond to that of a midwing position. Then use this value 

of Си to estimate a value of C;, that will apply to any wing 

position. 


Vertical tail.-— The contribution of the vertical tail to _ 


C, is usually estimated by the formula 


TOOL 


where Cru is preferably obtained from force-test data. 


As pointed out in reference 88 the esti- _ 
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Ftuvnx 183.—Charts and formula for estimating Cr, for the case of subsonic incompressible flow. Taken from reference 88. 
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When experimental data on Crp are available, the follow- 


ing formula from reference 88 can be used and will probably 
he more reliable than equation (30) because it takes into 
account any interference effects that might cause the effec- 
tive vertical location of the center of pressure of the tail to 
be different from the location determined by geometrical 


procedures: 
„ме. (с E il (31) 


or with the effective tail length —C,, ый йшй substituted 
for the geometric tail length l/b, 
дА 
т fiait Р 
1 Gee (=) LM (318) 


с, (25) oe 
Cr Цету ns 


C 


ьн 
1 


^ 

Wing.—The theory of reference 25 gives values of Ше 
derivative Cy (lateral force due to yawing) for the wing for 
а taper ratio of 1.0. The experimental data of references 25 
and 60 indicate that this theory is inadequate for making 
reliable estimates of Cy, It is recommended therefore 


that the experimental data given in references 25, 59, 60, 
and 61 be used in making estimates of Cy 


Vertical tail—The value of Cy, 


the formula 


Tein a 
can be estimated by 


А 
Ста РЕА b C Y Bratt (3 2) 


г 


or by Ше formula, in which the effective teil length 
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ui m / Cris is substituted for the geometric tail length #/6, 


C 290 
Y агр Ate it 


The discussion of lag-of-sidewash effects for Cs, | 
also to Су 


tail 
EFFECTS OF MACH NUMBER 


The effects of Mach number on the lateral stability deriv- 
atives have been treated theoretically in many investigations 
(see table IT) but very little experimental data have been 
obtained to verify this theoretical work. Moreover, only 
a small part of this experimental work has been covered in 
published reports (reference 114) because most of it is clas- 
sified at the present time. It appears, therefore, that 
estimates of the lateral-stability derivatives for the time 
being will have to be based largely on theoretical work. 

The effects of Mach number on the stability derivatives 
van be usually considered negligible for all airplane com- 
ponents except the wing and vertical tail. For the low-lift- 
coefficient condition in the case of many high-speed airplanes, 
the vertical tail contributes more than the wing to.all the 
stability derivatives except C; . For this reason, in ealcu- 
lations for transonic or supersonie speed conditions it is 
especially important to know the effects of Mach number on 
the vertical-tail lift-curve slope or Cr, ec 


Wing.—The effects of compressibility on the subsonic 
stability derivatives of the wing can be estimated by the 
formulas of reference 26. The values of the supersonic 
stability derivatives for some wing plan forms can be 
estimated by the references tabulated in table П. In this 
table the derivatives are grouped according to the type 
of wing plan form and to the particular derivatives covered. 
А helpful summary and discussion of the effects of Mach 
number on the derivatives for several different wing plan 
forms are presented in, reference 106. A summary of the 
theoretical lift-curve slope, damping i in roll, and center-of- 
pressure characteristics of various wing plan forms is pre- 
sented in reference 110. In the cases in which the theory 
shows large or abrupt changes in a stability derivative with 


apply 


ehanges in Mach number (for example, fig. 10 of reference. 


106} special care should be taken in estimating the deriva- 
tive in that particular Mach number range. The abrupt 
changes should be smoothed or faired out in a manner similar 
to that suggested in the following section for estimating 
Oy, а 

In some cases, experimental data for supersonic speeds will 
be available on the sideslip derivatives and on the damping- 
in-roll derivative C; . In such cases the experimental data 
should be used in preference to the theory. Some experi- 
mental results Бате indieated that the effect of the vertical 
location of the wing on Ше fuselage on the derivative Сү, 
might be greatly different at supersonic speeds from that at 
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subsonic speeds. Since no methods are presently available 
for estimating this effect for the supersonic case, it appears 
that, at least in the case of high-wing and low-wing designs, 
force teat data are necessary for obtaining an accurate 
estimate of Си. | 

Vertical tail.— The sideslip derivatives produced by the 
vertical teil at transonic and supersonic speeds can be esti- 
mated theoretically but should be obtained from force-test _ 
data whenever possible. These sideslip derivatives can be 
used to estimate the tail contributions to the other deriva- 
tives as pointed out previously. In estimates of the value 
of Сү, for transonic and supersonic speeds, corrections 
must be made for the effect of Mach number on the lift- 
curve slope of the tail, and these corrections should account 
for any differences in the end-plate effect of the horizontal 
tail on the vertical tail. 

For Mach numbers below about 0.8 or 0.9 and above _ 
about 1.6 or 1.8 the effect of Mach number on the lift-curve - 
slope of the vertical їа can be estimated satisfactorily from 
the theoretical values of references 26, 34, and 110. Since 
experimental data indicate that theoretical values of lift- ~ 
curve slope are usually too high for Mach numbers from 
&bout 0.8 or 0.9 to about 1.6 or 1.8, the empirically deter- 
mined fairings shown in figure 14 are recommended for use 
8S & guide in the use of the theory to obtain approximate 
estimates in this Mach number range when force-test data 
are not available. 

Experimental data have indicated that for vertical-tail 
configurations which have a tail length (distance from the 
center of gravity to the tail center of pressure) that is rela- 
tively short in terms of tail chords, the rearward shift of the 


tail center of pressure at supersonic speeds сап cause an  —— 


appreciable increase in the tail length and consequently an 
appreciable increase in the magnitude of some of the tail 


derivatives. Theoretical center-of-pressure positions for var- 


ious plan forms at supersonic speeds are given in reference 


110. : 
EFFECTS OF POWER 


On the basis of existing information, the effects of power 
on the lateral stability derivatives appear to be negligible 
in the ease of jet-propelled airplanes but these effects are 
often very large in the case of single-engine propeller-driven 
airplanes. Methods are available for estimating some of 


these power effects but in most cases experimental data are = 


necessary for making а satisfactory estimate. The effects 
of power can be broken down into two general classes: 
(1) The effects of the lateral force produced by the рго- 
peller itself | 
(2) The effects of the propeller slipstream оп the wing, 
fuselage, and vertical tail of the airplane 
Effects of propeller lateral force.—4À method of estimating 
the propeller-lateral-force derivative Cy, is presented in ref- 


егелсе 31 which is based on the work of references 32 and 33. 


= 2 o. 
Se -= F- 





FIGURE 14.—Examples of suggested fairing of theoretical values of lift-curve slope for use in 
estimating values for the vertical tall In the transonic range. 


The contribution of the propeller lateral force to the other 
stability derivatives can be estimated from this derivative 
by assuming that the ргорсПег is effectively а vertical tail 
surface and by using the expressions for the tail contribu- 
tion to the various derivatives presented in the preceding 
sections. Some experimental data on the effect of wind- 
milling propeller on all of the derivatives аге presented 4n 
reference 66. 

Effects of propeller slipstream.— The effects of propeller 
slipstream on the lateral-stability derivatives are usually 
much greater than the effects of propeller lateral force in 
the case of single-engine tractor airplanes. The slipstream 
effects on the wing, the fuselage, and the vertical tail can be 
considered as three independent effects. 

The slipstream effects on the wing can usually be neg- 
lected except for the derivatives C and C,. Experimental 
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data showing the decrease in effective dihedral (- Си) 
with power for single-engine airplanes are presented in ref- 
erences 55, 56, 57, 76, and 82. It appears highly desirable 
to determine this effect of power experimentally because 
interference effects make accurate estimations of the effect 
very difficult. The effect of the slipstream on the value of 
Ci, ‘cannot be estimated from the data on Ch 05 


described in the section on C. In fact, this procedure 
would probably give the wrong sign for the ineremenl of 
C, us - contributed by the slipstream. An approximation of 


this 1 increment might be obtained by estimating the slip- 
stream velocity and the lateral displacement of the slip- 
stream caused by yawing. Usually the power effects on 
C ы &nd Uus wil be greatest for the flap-extended 
configuration. 

In the сазе of the single-engine airplane the effect of the 
slipstream on the fuselage is usually to increase negatively 
the values of Cas and Cy,. (See references 55, 56, 57, 73, 
76, and 78.) Since no accurate methods of estimating these 
slipstr eam effects on C,, and Cy, are available, it is necessary 
to determine them from force-test data. 

The effects of the slipstream on the vertical tail are often 
very important and should also be determined from experi- 
mental data, if possible. The increase in dynamic pressure . 
at the tail caused by the slipstream is treated theoretically 
in reference 119 and is illustrated by the experimental data 
of references 51, 55, 56, 57, 73, 76, and 78. Tbe experimental 
data of referancé 78 also show thut the propeller slipstream 
can cause a destabilizing sidewash at the tail which will tend 
to reduce the stabilizing effect of the increased dynamic 
pressure at the tail. Since these date indieate that slip- 


stream effects оп the vertical tail vary greatly with airplane 


tion), use of experimental data appears to be the only satis- 


-factory estimation procedure at present. 


Suggested estimation procedure for power effects.—Thie 


following procedure is suggested for estimating power effects. 


Use tail-on 
Estimate rolling and 


Obtain force-test data for tail off and tail on. 
data directly for Сур Су, and Cry 
yawing derivatives as follows: 


(1) Estimate Cy, | from reference 31 and use this 
prepeccer 


derivative and proper linear dimensions to estimate other 
propeller derivatives (rolling апа yawing derivatives) in the 
same manner as tail derivatives. 

(2) Subtract tail-on data from tail-off data to get values 
of Cy, m ‚ Ca , and C, ra for the power-on condition and 
use ilicis AME to estimate the tail contribution to the 
other derivatives. 

(3) For tail-off values of rolling and yawing derivatives, 
use same values аз for power-off for all derivatives except 
C,. Estimate Са as suggested in preceding section. 

(4) Add the values obtained in steps 1, 2, and 3 to get the 
rolling and yawing derivatives for the complete airplane. 


INADEQUACIES IN PRESENT INFORMATION AND METHODS 


In the course of summarizing the estimation methods for 
the various stability derivatives, the need for much additional 
information on all the derivatives became apparent. In 
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particular, information is needed to aid in the estimation 
of the derivatives in the transonic and supersonie speed 
ranges. Additional work also needs to be done in correlating 
and analyzing existing subsonic data and in obtaining new 
experimental data for the development of semiempirical 
methods of estimating the subsonic derivatives without 
resort to force-test data. Another important need is for 
full-scale experimental results at all speeds for checking both 
low-seale data and the existing methods of estimating deriva- 
tives. Details of the need for additional work along these 
lines are discussed in the following sections. Studies should 
also be made to determine the conditions for which the use 
of steady-state stability derivatives in conventional stability 
equations is inadequate and to determine satisfactory 
methods of treating such conditions. 


TRANSONIC AND SUPERSONIC SPEEDS 


Additional theoretical work is needed on the estimation of 


stability derivatives in the transonic and supersonic speed 


ranges to cover the range of wing plan forms for all the 
derivatives. In particular, more work is needed on plan 
forms currently under consideration, such as wings having 
moderate sweepback and taper. This need is illustrated by 
table II which indicates that very little material is available 
on the stability derivatives for such plan forms except, 
perhaps, for the derivative C, . It appears from the table 


that this derivative and the triangular plan form have, in 
the past, received а disproportionate share of attention, 
probably beeause of the greater ease with which they could 
be treated theoretically. 

The greatest need for. work on stability derivatives at the 
present time is probably in the measurement of Ше deriva- 
tives at transonic and supersonic speeds. Experimental 
data on wings are urgently needed for checking the theoretical 
work and for use in the development of empirical corrections 
to the theory wherever necessary. Such corrections are 
particularly needed for fairing out abrupt variations of the 
derivatives with Mach number and for fairing through the 
Mach number range for which theory predicts infinite 
values. Examples of such discontinuities as indicated by 
theory are shown in figures 8 to 13 of reference 106. Since 
experimental data obtained at supersonic speeds on wing- 
fuselage combinations and on complete models have revealed 
interference effects that are different from those obtained &t 
subsonic speeds, it appears highly desirable to obtain at 
least a limited amount of experimental data at transonic 
and supersonic speeds to evaluate these interference effects. 
For example, investigations should be undertaken to deter- 
mine the effect of wing-fuselage interference on the derivative 
Си and the end-plate effect of the horizontal tail on the 


lift-curve slope of the vertical tail. 

Most of the experimental data on stability derivatives at 
transonic and supersonic speeds will of necessity be obtained 
at Reynolds numbers considerably less than full-scale values 
and under test conditions which might render the results 
open to question in some cases. Full-scale checks in flight 
of the low-scale data and of the estimation methods therefore 
appear to be desirable. Consequently the methods of meas- 
uring stability derivatives in flight now bemg developed Бу 


the Cornell Aeronautical Laboratory, Ше Massachusetís 
Institute of Technology, and the NACA should be extended | 
to transonic and supersonic speeds when the methods appear 
to be developed to a satisfactory degree of reliability for the 
subsonic case. Some preliminary considerations involved 


in the use of these flight techniques are discussed in, references | | 


120 to 123. 
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SUBSONIC SPEEDS eMe 4o 


The methods presented in this report for estimating the 
stability derivatives at subsonic speeds depend either directly 
or indirectly on the use of force-test data. These methods 
are probably more reliable than methods which do nob _ 
involve the use of force-test data on the particular design - 
under consideration ог on а similar design. Methods which 
do not rely on such data are desirable in some cases, however, 
because the necessary data wil not always be available. 


In the case of sideslip derivatives, empirical methods can _ 


probably be developed largely from existing information. 
In some cases ii wil be necessary to augment the existing 
information with new results since much of the available 
force-test data were not obtained in a manner that would 
make the data readily usable for developing general esti- 
mation procedures. 

In the case of rolling and yawing derivatives, considerably 
jess information is available than in the case of the sideslip 
derivatives. Most of the information now available was 
obtained in the Langley stability tunnel, principally on wing 
configurations and to a limited extent on complete airplane 
models and airplane components other than the wing. 


Considerably more work js required, especially for com- 


ponents in combination, before satisfactory methods can be 
developed for estimating rolling and yawing derivatives 
without the use of force-test data on the particular design 
under consideration or on a similar design. г. 

In discussing the work necessary for developing new рго- 
cedures for estimating the stability derivatives without the | 


use of force-test data on the design under consideration ог 


on а similar design, it is useful to break the problem down > 


--- 


into two parts: (1) effect of individual components and (2)  — 


the effect of interference of the components on each other. 


The principal components to be considered are the fuse- — 


lage, wing, vertical tail, and propeller. For the isolated 
fuselage, the main problem is the development of methods 


for the estimation of Са, апа then, perhaps, of C,, and 
Cy, For the isolated wing, the main problem is to estimate 
the derivatives at lift coefficients above that at which 


separation begins. Such estimations сап be made with 


reasonable accuracy for some of the derivatives by existing | 


methods which make use of force-test data, but the develop- 
ment of methods which do not involve the use of force-test 
data will probably be very dificult. For the isolated vertical 
tail, the problem is to establish the effective teil area and 
aspect ratio from the geometry of the tail so that the lift- 
curve slope (or Су,) of the tail can be caleulated. Solutions 
to this seemingly simple problem have in the past become 
involved with interference effects so that, аз yet, no reliable 
methods have been published for estimating Су, of the 


vertical tail from its geometry. For the isolated propellers, 


850 


the work that is needed at present is & systematic check of 


existing methods of estimating the lateral force on the 


propeller to determine the accuracy of these methods. 

The principal interference effects to be considered are 
mutual interference of the wing and fuselage, wing-fuselege 
interference on the vertical tail, horizontal-tuil interference 
on the vertical tail, and propeller-slipstream interference on 
the wing, fuselage, and vertical tail. The mutual-interference 
effects of the wing and fuselage are probably important 
only for the derivatives Cu C," апа С,. А large amount 
of experimental data is available for the sideslip deriv- 
atives but no procedures for estimating the interference 
effects on these derivatives have been reported. Wing- 
fuselage interfefence has very important effects оп Cy, of 
the vertical tail and consequently on all of the stability 
derivatives for some flight conditions. These effects result 
from the sidewash and change in dynamic pressure at the 
tail which may result from sideslipping, rolling, or yawing. 
Although considerable data which show these interference 
effects are available, particularly for the case of sideslipping, 
no reliable methods exist for estimating the interference 
effects. Horizontal-tail interference also has an important 
effect on Cy, of the vertical tail for some horizontal-tail 
positions. Some work on a limited number of configurations 
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has been done toward developing methods of estimating 
this effect. but data are required on more configurations 
before the generally applicable methods can be evolved. 
The propeller slipstream can cause important effects on 
C ‘and C, of the wing, on Ca, and Cy, of the fuselage, and 
on Cy, of the tail (and consequently on the tail contribution 
to all the derivatives). Some data are available for the 
affect of the slipstream on the sideslip derivatives but, 
because of the complexity of this problem, considerable 
additional data may be required before a satisfactory method 
of estimating the slipstream effects can be developed. 

As mentioned in the preceding section, full-scale checks 
of low-scale data and of the estimation methods are desirable. 
For the subsonic case some of the checks can be obtained 
from large-scale wind-tunnel tests but some cheeks in full- 
scale flight tests should also be obtained when the various 
methods of measuring stability derivatives in flight have 
been developed to а satisfactory degree of accuracy. 


LANGLEY AERONAUTICAL LABORATORY, 
NATIONAL ÁpvisORY COMMITTEE FOR AERONAUTICS, 
LANGLEY FrELD, Va., December 18, 1950. 


APPENDIX А 


EQUATIONS OF 


LATERAL MOTION 


The dimensional equations for the lateral motions of an airplane are 


Фф dL do dy ody dL y 

mes" df Әр dei m mk xs ДЕ d& or dt Ov” —L.=0 (A1) 
d aN d$, dy Ody ON мо. 

mExz dà op di "tx d? or di Ov М„=0 (А2) 

ms SP (if) фт VR e - 2 р. —, —(Lift)(tan y) --m 2 ОГ. v—Y.=0 (A3) 


If equations (A1) and (A2) are divided byz p V*Sb and equation (АЗ) is divided by = p V*S, the equations of motion may be 
expressed in the conventional nondimensionel form in which they have generally been presented in NACA papers (for 


example, see reference 2): 


1o dé a? 
wks rfr C Get Dua qi 
ык 10, Ben 

1 de 1 


-а Су 19 —Cró--2u 


2 р 


1 
2 
1 
2 


© 


10, —О„8—Со„=0 
dy 
s. ds Cw —ÜC,.— 0 (А 4) 


5 Ст, ФУ. Cz (tan Y)v--2u Borg- Cy, =0 


In order to convert these equations into a form which will reduce the number of arithmetical and algebraic steps in 
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performing stability calculations, equations (А4) are multi- 





„Каз K 














: : : K K—— 
plied by m/pSb and are written in the following form: d Ку? EK? ` 
(Г Рф-(АЕ-ТИ) у — 188 —l.=0 L=- C, па= Е б eh Cy 

a 2-и 
(K,D'—n,D)6-4-(D* —n,D) пав —n,—0 2Ky 2K; в 2 a 
Cr С. GS жы Б muse Y 
(750-52) $4(D—y,D—— ten r) + "OAK; Y Кр в Уат 
(D—y4) 8—y.—0 — =l A 
where l; 4Ky Ci, Th 4K, Cs. Yr m Cy, 
em Lom ut _ 4 __А и і 
„== бъ PS e D—q.| SIRAC ое9куб. Yge 
APPENDIX B 


APPLICATION OF THE LAPLACE TRANSFORMATION TO CALCULATING MOTIONS 


The application of the Laplace transformation to the 
calculation of the lateral motions of airplanes is presented in 
order to illustrate the development of the equations of motion 
in the form in which they are presented in the present report. 
This work is similar to that presented in references 5 and 6. 
In fact, it follows the presentation in reference 5 very closely. 
Reference 6 presents a brief explanation of the Laplace 
transformation and its application to solution of the equations 
of motion of an airplane. This report also makes reference 
to detailed explanations of the Laplace transformation. In 
cases where modification of the equations presented in the 
present report are necessary, reference should be made to 
these texts for an understanding of the mathematics involved. 
Applying the Laplace transforms 


L=} Цр 


L($)— dx L(D*6)— Md, —Xóy—(D д} 


&nd multiplying each of the equations by X transforms 
equ&tions (АБ) from appendix А to 


(МФ (КА — LX, ВАВ 

(Ka —п„А)ф, (№ — п, Ay — ng By — r3 

сета дере баа y» | vat 
(A — ув\)В›==т» 


(B1) 


where 


пя „МФ +(К.№ 1А) MD ФУ KMD yjet l: 


г (Ко —n5X)dot- A — n Мо KAD pjt A Dpt ne 
з-- —Yprdo +A—Yy A) por ЖВа-- Ye 


Solving equations (B1) by determinants gives 


— lg Fi KX -LN 

— fig n A3—n,A* 

Мул от toy SZ (tan y). 
маи ЕА 1А? МТА? 

—ng KA —n,M МАЗ 

A*— yar уум 62 Х Mya СЕ (tan у)А 
which may be expressed as 


а а -- а + (4A? -- аА +a, 


h=- Ant Be Ot DLA 82) 
Similarly, the expressions for уу and В, are 
pa БА E BEBO EDO E BAF BF (B3) 
С AXAMLBN+AON+ DALE) 
ER Cot H е,\%-Е COA HGA -- C, (B4) 
Bi МАМ + B3 H-CX*4- DATE) 


where the expressions for the coefficients in equations (B2) 
to (B4) are given in terms of the mass and aerodynamic 
stability derivatives by equations (1) to (4) in the main body 
of this report. 

In order to obtain the actual variables from the trans- 
formed variables, an inverse Laplace transformation must be 
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applied. The expressions for фу, Ya, and f, are the form 
ald where ux and в, are polynomials, the degree of » being 
higher than that of v. The inverse transform of a function 


of this type is 
L^ о? ЩА). on, 
nml 


бу 2 (А) Б 5) 





In this equation all of the roots А of m=0 are assumed to 
be distinct. This assumption is valid for fa; but for фу and 
vy, 2=0 has two zero roots. (See equations (B2), (B3), 


and (B4). The terms in the equations for ¢ and y resulting 
from the two zero roots are 
dQ 
= (0)--9(0)« (B6) 
where 
Q=) 
РА 


From equations (B5) and (B6) then, Ше inverse transforms 


of d», у, and 8, (equations (B2), (B3), and (B4)) are 


ф-- d, HAt дем 4 A e H- Дье Aa (B7) 
' y= Вем+ Bitit Bst Ве M-- Во + В. (B8) 
B= Сем+ОСемз+ Сит Сем - C, (B9) 


The equations for the rolling velocity p and the yawing 
velocity r can be obtained from equations (B7) and (D$) Бу 
differentiation 


p=ž (Ame HA MF Ae HA ЧА)) (B10) 


r= (Biens + Be ™ 4- Вме”з+- Вмеч+ В.) (Ві 1) 


where the expressions for the coefficients of equations (87) 
to (B11) are given by equations (6) to (8) in the section 
entitled “Calculation of Motions.” 


APPENDIX С 
SOLUTION OF BIQUADRATIC EQUATION 


Many methods are available, of course, for solving for the 
roots of в biquadratic equation. For example, there are 


Horner’s, Ferrari’s, Bernoulli’s, Descartes’, and Hitcheock’s | 


methods; various methods of solution by trial; and also 
various graphical methods such as that given in reference 1. 
Solution by trial in which synthetic division is used, however, 
is recommended as being the simplest method for most 
lateral stability work. The characteristic equation for the 
lateral motions of an airplane 


АМУ BMN+C¥+DA+E=0 


generally has two real roots and a pair of conjugate complex 
roots. For these cases the two real roots can be factored out 


easily and the remaining quadratic solved for the conjugate . 
complex roots. In the few cases for which all four of the. 


roots of the characteristic equation are complex, Descartes’ 
method can be used to factor the biquadratic equation into 
two quadratics. When there аге real roots, solution by 
Descartes’ method requires more time than factoring out 
the real roots singly and consequently is not recommended 
for general use. These methods of solution are explained 
in the following sections. 


SOLUTION BY TRIAL BY MEANS OF SYNTHETIC DIVISION 


Solution for real roots by trial by means of synthetic 
division consists of successive approximations of а root and 
checking by synthetic division until the root is determined 
to the desired degree of accuracy. This check by synthetic 
division is based on the fact that if a is а root of а polynomial 
F(z) then z—a is a factor of f(x) and consequently no re- 
mainder is left when f(x) is divided. by z—a. 

The method of solving the stability biquadratic equation 
by trial with synthetic division is explained in three steps in 
the following sections. First, the rule for synthetic division 
and a numerical example are given. Second, the specific 
use of synthetic division for factoring a biquadratice is illus- 
trated by a simplified example for which the roots are known. 
This example shows how the cubic and then the quadratic 
factors of the biquadratie are obtained. "Third, the use of 
synthetic division in extracting the roots of a representative 
characteristic stability biquadratic is illustrated with special 
reference to methods of making the first approximations of 
the real roots. | 

Explanation of synthetic division.— Synthetic division is 
explained in almost ali algebra text books but is presented 
herein for the convenience of the reader. The rule for 
synthetic division may be given as follows: 

Assume that a polynomial in г (f(z)) is to be divided by 
—а; write the coefficients of the polynomial in order, 
supplying 0 when а coefficient ts lacking. 


Multiply a by the first coefficient, and add (algebraically) 
the product to the next coefficient. 

Multiply this sum by a, add to the next coefficient, and 
proceed until all the coefficients.are used. The last sum is 
the remainder and also the value of the polynomial when a 
IS substituted for the variable г. 

For example, divide z*-- 323 -3:?—2z—6 by z—3 


1+3+ 3— 1— 6 
+318463 +186 | 3 | 


Use of synthetic division in factoring out roots.—The use 


of synthetic division to factor out two known rational roots 


of а biquadratie equation. is illustrated by the following 
simple example. 
real roots of the characteristic stability equation which, of 
course, are not normally known but can be approximated by 
the method given in the next section of this report. 


One factor of the biquadratic is x—1 so there is no re  — 


mainder when the biquadratic is divided by the root 1 
1+3+3—i—6 


+1+4+-7+6]1 
1+4+7+6 0 


Since the remainder is 0, x—1 is one factor of Ше biquadratic 
equation and 2*?-++427+72z+6 is another factor. 
as а cubic equation must have at least one real root, a second 
real root of the biquadratic equation can be factored out 
of the cubic. For example r+2 is a factor so divide the 
cubic by the root —2. 


“дал oO 


1+4+7+6 
—2—4—86|—2 
1+2+3 0 


The factors of the biquadratic then аге 2—1, x+2, and 
z1--2z-38. 


roots by the quadratic formula. For example 


-— l4 .19 
r= Tivi LERZ 


Example of application to characteristic equation.—Rea- 
sonably accurate first approximations to the real roots of the. 
characteristic equation can be obtained from simple formulas, 
Successively closer approximations can then be obtained by- 
interpolating from the remainders. The following ex&mplé 
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These two rational roots represent the two _ 


Inasmuch 


The quadratic factor can be solved for its — 
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illustrates the application of this method to obtaining the 
roots of the stability biquadratic. "The biquadratic 


M- 10.43? -1- 16.324? 4- 68.04 — 9.10 —0 


is of the form 
AM 4- B+ CX --DX4-E—0 


Since the coefficient Е is generally much smaller than 
coefficient D in lataral stability work, one of the real roots 
(usually the smaller of the two) is approximately equal to 
-ЕП or it may be more closely approximated by the 
equation 


| Е | 
те 08 
D 
or for the particular case 
—9.10 "T 
M aso 632—910) 029. 
` 68.6 


Approximating the root by synthetic division 


1+10.43+16.32+68.6—9.t0 Approximation 
+ .184- 1.304- 2.3+9.10|.1284 2 

+ 134 1.364 2.34+9.14|.129 1. 
14-10.564-17.68470.94- .04 1 
14-10.56 4- 17.68 4- 70.9 4- 0 | 2 


For this root, the second approximation was determined by 
dividing the coefficient E by the fourth sum from the quotient 


| —9.10 
70.9 


This procedure generally provides а good second approxi- 
mation for the small real root. 
The cubic equation obtained by setting 


+-10.564?-+17.68A+ 70.9 
equal to zero is of the form 
аз b-+ch+d=0 


In most lateral-stability work, a real root of this equation 
will be approximately equal to —Ё or it may be more closely 
approximated by the equation 


or for the particular case 


(10.56) 4- 70.9 


—(10.5693-17.68 9-6? 


у= 
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Approximating the root by synthetic division 










14-10.56 +17.68+70.9 Approximation 
— 9.48 —10.20— 70.9 | — 9.485 6 
— 9.49 —10.16—71.4 | — 9.49 5 
— 9.48 —10.25—70.4/-9.48 С 4 
—-9.45 —10.50—07.9 | — 9.45 3 
— 9.55 — 9.64— 76.8 | —9.55 2 

| — 9.65 — 8./8—85.9|—9.05 — 1 — - 
1+ 0.91 + 8.00--15.0 1 

i+ 1.01 + 8.04— 5.9 | 2 
13-141 d 4183-30 3 

‘1+ 1.08 + 7.48+ 0.5 4 

I+ 1.07 + 7.52— 0.5 5 

1+ 1.0754 7.48 0 6 


For this lerge real root there is no simple method of deter- 
mining the second approximation as there was in the case of 
the smaller real root. The magnitude of the estimated root 
in this case is arbitrarily increased or decreased slighUy from 
the first approximation. From the remainders determined 
from the first two approximations, a fairly closo third ap- 
proximation can then be made. 
Factoring the quadratic equation obtained by setting 


N+1.075A+-7.48 


equal to zero by use of the quadratic formula gives the final 
two roots of the biquadratic equation 


1.075 +4/1.16—29.92 
ee ха 


-— m - 


--овзва „86 | + 


= — 0.5384 2.682 


The roots of the biquadratic equation may be checked by 
multiplying the four factors to determine whether their 
product equals the original biquadratic 


(4 —0.1284) (1+ 9.485) (A+0.538 + 2.682) (A+ 0.538 — 2.681) = 
M+ 10.4308-+ 16.320? 68.6A—9.10 
SOLUTION BY DESCARTES’ METHOD 


Descartes’ method of solving a biquadratic equation is 
particularly useful for solving equations which do not have 
any real roots. This method is explained in most text 
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hooks on advanced algebra and theory of equations. In 
general, the method consists of reducing the biquadratic 
equation to а cubic equation which can be solved easily. 
One root of the cubic equation is used to form two quadratic 
equations the roots of which are used to obtain the roots of 
the biquadratic equation. 

Method.—Reduce the general biquadratic equation 


895 
end form the equation 


1 
+5 eis qi— т фрита 0 


and solve this cubic equation in zë for one of Из roots 2240. 
Solution by trial by means of synthetic division is recom- 
mended. Determine the values of l and m from the equa- 
tions. 


AM+ B+ CX EDAX4-E-0 (=i — ща 

to the form T" id 2 ir 

M+ 53-4 - eM дА Ее-- 0 -9 Ent. 

тота T 
by dividing by А. Substitute the values of / and m and the value of x used in. 

Obtain the values of q, r, and 8 from the following equa- obtaining Г and m in the equations u 
tions: до 
I y +2ry+i=0 | эс и 

а==б—с р? y?—2ry-+m=0 | 


and solve these quadratic equations for their roots y from 
which the roots of the biquadratic equation may be obtained 
from the following relation: 


duel b? 


г. T 3 


b 
4 — ра 
F16 258 E 
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APPENDIX D 
SPECIAL NOTATION USED IN CALCULATING MOTIONS WHEN THE CHARACTERISTIC EQUATION HAS COMPLEX ROOTS 


When two of the roots Х and X; are conjugate complex, the coefficients А, and Az B; and B}, C, and Ca will be conjugate 
complex. If R--IH is one of the roots ^, and if the powers of А, are expressed as 


—R,4 М 
=A iy 


¥ 


then 


AP = Rat Lot 
№= Ryt dh 
=R;+ Iċ г... 
ЕА: 
Substitution of the root А-Л in the expression for .4; gives 


(aoB + а а, Ез-- аз F3 - ай: + a4) t (a7, - a1143- а-а: Р, tahi 
Д5, 1-58. 440R +3 DR: 42ER) (ВАТ, -5ВІ, +40 h+3 DhE) 


The division of these complex numbers is indicated by the equation 


21171 


ог Буше, Yayi— 2102. 
2122001021 TYayi Xia 
га узо T 


г |-у» ха уз? 


It в evident from these relations that A; is а complex number. In this case new symbols аге used to represent the геа] арда — 
imaginary parts of А, as follows: 
я Аз RAT IA 
‚4, is the conjugate of 4, and will be referred to as 
A= Ra— Га? ` 
By procedures similar to those for the A coefficients, 


p, (os Bak bus bRat bit bt ВЕБ LE bat baht БАТ 
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which may be referred to as 


aud 
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В„=Ё„— 


Also, 


(Coat cR аз На Вос В1) Cols сев teala telji 


Веви. 


Isi 


~ GAR + SBR. + ACRES DB, V 2ER)--(6À 1,4- 5BI,4- 4C I,2-3DI T 2E Li 


which may be referred to as 


and 


Similar analysis shows that, if the roots M and M are also conjugate complex quantities (X,— R'--I^i and X, ^— 


and 
A= R' АТ. 
where 
re («Вай Без Ва е Е .-- aR’ tH atla" а Haal’ tH ata'i TN 
(GAR. +5 R..4+4CR',+-3DR.42ER')+(6Al' + 5BI',.--4CI',--8DI' 3-2 EI'yi 
Also, . 
== Ff’, +I’ pi 
and 
By= R'y—I' м 
where 
(boR s+ 0.2", 6,0% РОСА И иа БВ bal ЕВА ВВ) 0 — 
С (6AR’,+5BR .+4CR,4+-3DR+2ER)+(6Al + БВГ A402, +ЗРГ.Ф2ЕГ т, 
Similarly, 
C,—R'ec4- I'cà 
and 


C = Re + Fei 


C= Во — 


Ici 


ft), then 


Ay ВГА 









C,— R'c—I'ci 


where * 


4. 


~t 


C= (cot, - c B E eR at eR eR Hle 


(GABY 5BR',--4C 
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TABLE I.—TABLE FOR CALCULATING OSCILLATORY STABILITY BOUNDARIES 


CONSTANTS 


b, feet 
S, square feot 


= (7)-(8), deg 
Kx! =[(10)/(1)] 1сот(18+-(П)/Суряги (12) 


| 4Кх!=4х (13) 


Kg! -[(1:2/0)]1609102) -[(10)/(0 Tri n1(12) 
«Ку e AX (18) 
Кхк=(1(11)/0В—[(10)/(1)})ай(12) eos (12) 
Кл (17/019) 

Къ (7/16) 

ЮК:-- (18) X (19) 


: me (8)/182.2Х(1) X (2) X (5)] 


du 4X (21) 


. (Ci/2) tan += (0) X tan (0) 
А =1-— (20) 
: Ka (СД) tan у= (28) X (19) 


Ka (C12) tan у- CA (35) --(6) 
2x um 2X (18)/(21) 



















Total derivatives 


(49) 


COMPUTATION PROCEDURE 


. Fill in values for the known geometric and mass 


characteristics and flight, conditions (constante 


(1) to (11)). 


. Determine v&lues of the tail-off stability deriva- 


tives (constants (28) to (35)) by methods pre- 
sented in the text. 


. Select &nd enter in column (36) values of the inde- 
pendent уапа е C,, to cover the range for 
which the boundary 1s required. 


. Work out columns (37) &nd (38) to obtain velis 
of Суда" 


. From these values of Cr, ca? determine the tail 


contribution to the stability derivatives (col- 
umns (39) to (44)) by methods presented in the 
text. 


. Perform the operations indicated for constants 


(12) to (27) and columns (45) to (90) to obtain 
the values of Routh’s discriminant R. 


. Solve for values of и the quadratics formed by 
setting the values in column (90) equal to zero ` 


(columns (91) and (92)). 


. Solution of columns (93) and (94) gives the values: 


of Са, required for neutral oscillatory stability 
for the values of C,, in column (36). 





Oombined mass and aerodynamic derivatives 
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(88) (82) (85) (80) 
CL Ог 
С, jz tan y | Са g n y 
C — 2 Xin E 
„| ae Der б, 
(36 (ох) | (29) х(Б8) | (8x (66) | (28) (74) — (88) |, 
УР) Даде ки 
9 —, 0083 —,0110/a 
-+0 
„05 — 0471 - ДИ 
—, 8418 
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в 1 =, 1245 
А Burt 
» 25 —, 2090 —. 19632 
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(87) (88) 
вор AID 
(m)XCDX(8) — | сюжвохви 
02781 , 1380 
Na- 0010 А пич 
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АА ав, 0008 ИДИ А ад, 0508 
шүү Тт" 
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— 168. 9445 


Съ, И а ВТ. 113 


„1838 
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—, 10811 
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(93) (94) 
Си Ci, 
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. 2008 —, 0770 
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TABLE IL—REFERENCES CONTAINING USEFUL INFORMATION FOR ESTIMATING LATERAL STABILITY DERIVATIVES 
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SEEEN 
ар&®& 
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ВЕЕР 


~J = 






, 30, 40, 41, 42, 43, 44, 
45, 46, 47, 48, 49, 50, 
Bl, 52, 53, 54, 55, 66, 
57, 58, 59, 00, 61, 62. 
63, 64, 65, 80, 67 

88, 39, И 41, 42, 42, 44, 

, 48, 41, 48, 49, 50, Х 

51, 52, 53, 54, 55, 56, 
57, 58, 69, 60, 61, 62, 
63, 84, 65, 68, 67, 73, 
74, ТА, 76, 77, 78, 79, 
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98, 100 
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Supersonic (аП are theoretical estimation methods) 
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